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ABSTRACT 


To determine the feasibility of coupling the output of an optical fiber to a rib 
waveguide in a temperature environment ranging from 20°C to 300°C, a theoretical 
calculation of the coupling efficiency between the two was investigated. This is a 
significant problem which needs to be addressed to determine whether an integrated 
optic device can function in a harsh temperature environment. Because the behavior of 
the integrated-optic device is polarization sensitive, a polarization-preserving optic 
fiber, via its elliptical core, was used to couple light with a known polarization into the 
device. To couple light energy efficiently from an optical fiber into a channel 
waveguide, the design of both components should provide for well-matched electric 
field profiles. 

The rib waveguide analyzed was the light input channel of an integrated-optic 
pressure sensor. Due to the complex geometry of the rib waveguide, there is no 
analytical solution to the wave equation for the guided modes. Approximation or 
numerical techniques must be utilized to determine the propagation constants and field 
patterns of the guide. In this study, three solution methods were used to determine the 
field profiles of both the fiber and guide: the effective-index method (EIM), 
Marcatili's approximation, and a Fourier method. These methods were utilized 
independently to calculate the electric field profile of a rib channel waveguide and 
elliptical fiber at two temperatures, 20°C and 300°C. These temperatures were chosen 
to represent a nominal and a high temperature that the device would experience. 
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Using the electric field profile calculated from each method, the theoretical 
coupling efficiency between the single-mode optical fiber and rib waveguide was 
calculated using the overlap integral and results of the techniques compared. Initially, 
perfect alignment was assumed and the coupling efficiency calculated. Then, the 
coupling efficiency calculation was repeated for a range of transverse offsets at both 
temperatures. Results of the calculation indicate a high coupling efficiency can be 
achieved when the two components were properly aligned. The coupling efficiency 
was more sensitive to alignment offsets in the y direction than the x, due to the 
elliptical modal profile of both components. Changes in the coupling efficiency over 
temperature were found to be minimal. 
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CHAPTER I 


INTRODUCTION 

1.1 Fiber Optic and Integrated Optic Benefits 

Transmitting fields via lossless dielectric cylinders was theoretically demon- 
strated as early as 1910 [1]. This concept was not experimentally studied until the 
advent of low-loss optical fiber [2], which has revolutionized the communication indus- 
try. Copper telephone wires are now being replaced by optical fibers due to benefits 
which include: fewer repeaters required compared to electrical wires over the same 
distance; larger amounts of information transmitted over smaller diameter fibers due 
to the higher bandwidth; and an improved signal to noise ratio. 

The advantages seen in communications can be applied to aeronautics. On 
an airplane, optical sensors can be located remotely from the processing electronics 
with a fiber connection. Since optical sensors can be made compact, they can non- 
intrusively detect environmental parameters. These sensors are inherently immune 
to electromagnetic interference (EMI), such as radar, which can corrupt an electrical 
sensor’s performance. This is a very important property because sensors are often 
required to operate in harsh electromagnetic environments. To use electrical sensors 
in these environments requires additional components for proper shielding. For aero- 
nautics applications, this added shielding increases the weight of an aircraft, which 
is an undesired result as added weight reduces the amount of passengers and cargo 
that can be transported as well as the aircraft range and speed [3]. Unlike electric 
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wires and sensors, optical fibers and sensors do not exhibit sparking, short circuit, 
or fire hazards caused by the degradation of the wire insulation [4]. These hazards 
can present serious safety concerns on an aircraft. Ground loop problems are also 
avoided using optics which reduces noise. Transmitting sensitive information via op- 
tical fibers can improve security because the electromagnetic field strength outside a 
fiber is much smaller than that outside an electrical wire, thus making interception 
nearly impossible. Optical sensors are amenable to multiplexing, which can greatly 
improve efficiency. Sensors distributed in a matrix can map a measured physical 
parameter over the aircraft. Therefore, to replace electronic components with opti- 
cal components would have wide-spread benefit. However, because this technology 
is still maturing, it is more expensive than its electrical counterpart and hence not 
widely utilized. 

For aeronautics applications, an optical sensor must perform to its design 
specifications while located in a hostile environment. This is a challenge because if not 
properly designed, the high temperatures of an aircraft engine can seriously degrade 
or destroy the device. Modeling has become an important tool in predicting the 
waveguide’s behavior with temperature. To model the behavior of an optical device 
over temperature, changes in both the refractive index of the layers and the physical 
dimensions must be determined. Temperature-induced changes in the materials due 
to thermal expansion and stress must be taken into account in order to determine 
the device properties at an elevated temperature. 
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This research explored the coupling efficiency between a single-mode optical 
fiber and a rib waveguide. Single mode sensors are generally more sensitive than mul- 
timode sensors, which was one of the main reasons it was desirable to determine the 
coupling feasibility between single-mode components. The optical fiber transmits the 
source light into the guide, which serves as the input to an integrated optic pressure 
sensor. Due to the geometry of these two components, approximation and numerical 
techniques have been used to determine the electric field profiles and hence the theo- 
retical coupling efficiency. Because this sensor is to be used in a temperature varying 
environment on an airplane, changes in the coupling efficiency over temperature are 
of interest. 

1.2 Dissertation Objectives 

Presently, single-mode fibers and waveguides are not considered viable for 
use in an aircraft environment due to the coupling problem [5], [6]. The majority 
of devices are not optimized for coupling, resulting in reduced maximum coupling 
efficiency, such that a slight misalignment drastically reduces the power transfer. 
Thus, the results of fiber to rib guide coupling are a critical factor in predicting 
device feasibility because a sufficient amount of input light power is required for 
the waveguide output to be detected. If theoretical calculations indicate a high 
interconnect loss, then modifications must be made prior to fabrication with the goal 
to improve the coupling efficiency. 

Calculating the theoretical coupling efficiency between an elliptical-core op- 
tical fiber and a rib waveguide at two temperature conditions with various alignment 
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positions, to determine interconnect feasibility, was the objective of this research. To 
do so required the calculation of the electric field profiles of both components in order 
to evaluate the well-known overlap integral. Because there is no analytical solution 
to the field of the rib guide, two approximation and one numerical technique were 
used and the results of these techniques for the rib were compared. The two approx- 
imation techniques chosen, the EIM [7], [8], [9], [10], [11], [12] and Marcatili’s method 
[10], [13], [14], [15], are widely accepted techniques that have been extensively stud- 
ied. The numerical technique, the Fourier method, was chosen due to its reported 
excellent agreement with known exact solutions [16]. 

The fiber field can be calculated analytically using a series of Mathieu func- 
tions [2], however, the solution cannot be determined in closed form because the 
eigenvalue equation involves infinite determinants [17]. It was determined that the 
same methods used to calculate the rib guide characteristics were applicable to the 
fiber [16], [17], [18]. 

The effective index. Tie, related to the modal propagation constant, is a 
result of the field calculation and was used for comparison purposes. Both the 
and coupling efficiency calculations were performed at the two temperatures. As a 
consequence, it was necessary to determine the index and physical dimensions of both 
components at the two temperatures, producing a different field profile for each. 



5 


The two temperatures were chosen to represent a nominal temperature and 
a high temperature that such a device would experience near the inlet of an aircraft 
engine [4]. The maximum theoretical coupling efficiencies at each temperature were 
compared to determine if the temperature change was significant to affect the sensor 
performance. 

1.3 Organization 

The organization of this dissertation is briefly outlined here. Chapter II 
provides an overview of documented research related to this dissertation. A brief 
description of the sensor whose input is being studied is also given in Chapter II. 
The fiber and sensor physical parameters are discussed in Chapter III along with the 
physical parameters’ behavior with temperature. In Chapter IV the three methods 
used to calculate the electric field profiles of both the fiber and rib waveguide are 
discussed. Chapter V is divided into two parts, room temperature (20°C) and 300°C 
results. These results include the field profiles and corresponding coupling efficiency 
results for transverse offsets. Using the three methods independently, the results 
were compared. A discussion on maximizing the theoretical coupling efficiency by 
altering the waveguide parameters is given in Chapter V. The conclusion and plans for 
future work are discussed in Chapter VI. A list of acronyms used with corresponding 
definitions is provided in the on page xv. 



CHAPTER II 


BACKGROUND INFORMATION 


2.1 Introduction 

This chapter describes previous research and is divided into rib, fiber, and 
coupling efficiency subsections which provide an overview of the research regarding 
field solutions required to calculate fiber to guide coupling. The first two subsections 
explore the field profile calculations documented to date on the two coupling com- 
ponents. The third subsection discusses research on the coupling efficiency between 
a fiber and a rib waveguide. A few specific examples from the literature survey will 
be explored to illustrate the differences between previous research and that of this 
dissertation. A brief description of the pressure sensor for which the rib waveguide 
is the light input is also given. 

2.2 Previous Research 

The problem of transmitting the output light from a fiber into an integrated 
optic device is relatively new. Until the advent of low loss fibers in the 1960’s, it was 
impractical to attempt coupling, as the fiber exhibited such high loss. The commu- 
nications industry has addressed this coupling issue extensively for lithium niobate 
(LiNbOa) devices. Results from the literature search performed in preparation for 
this research will be discussed. 

The majority of the research calculates the propagation characteristics of ei- 
ther an optical fiber or a waveguide. For this calculation, there are many techniques 
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from which to choose. These methods include: finite element method, EIM, Mar- 
catili’s method, Fourier method, variational method, and beam propagation method 
(BPM). 

2.2.1 Rib Field Calculation 

Due to the complex geometry of the rib waveguide, approximation or nu- 
merical solution techniques are necessary to determine its propagation characteristics. 
Two review articles, one by Chiang [19] and the other by Saad [20], discuss the meth- 
ods that have been used to solve for the propagation characteristics of waveguides. 
A discussion of several of the methods will be given. 

The EIM has been used to calculate the profile of a rib [8], [11], [19]. Caution 
must be exercised when using this technique, because inaccurate results may occur if 
applied to a device which does not meet the approximation assumptions. For exam- 
ple, field boundary conditions are not satisfied using this approximation technique 
[10]; thus narrow ribs are not accurately modeled using this method. Determination 
of the accuracy of this method relies upon comparison with a numerical technique 
[21]. The EIM is known to be accurate for ribs with small etch depths, such as 
the rib tmder study. However, it consistently yields propagation constants that are 
larger than those provided by the numerical results. Because this technique is only 
applicable to a specific rib geometry, it is not widely recognized as useful for field 
calculations, but rather only for the determination of the propagation constant [10]. 
However, this remains one of the most popular techniques for the analysis of rectan- 
gular structures. Variations on this method are numerous [22], [23], [24], [25], [26] with 
the intent to improve results. 
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Marcatili’s method [13] has been applied to both rib and strip-loaded waveg- 
uides. Results from this technique compared well with numerical solutions [14], [15]. 
This method assumes the field is tightly confined to the core such that the field in 
the outer regions is negligible. The technique has even been applied to modes in 
anisotropic rectangular waveguides [27]. A perturbation correction to the propaga- 
tion constant calculated using Marcatili’s method can improve the accuracy [14]. A 
perturbation correction can also be applied to the EIM [28]. The major advantage 
of the EIM and Marcatili’s technique is their simplicity. 

Goell [29] reported a numerical technique based on the expansion of circu- 
lar harmonic functions. This method produces accurate results for the propagation 
constant, but becomes cumbersoipe for devices with aspect ratios other than unity. 
Many researchers refer to Goell’s work for comparison, yet it is infrequently used. 

The finite element method is an application of the Rayleigh-Ritz procedure 
to a set of basis functions which are confined to a small spatial region [30]. When 
applied properly, this method has been determined to be very accurate and can be 
applied to many geometries [31]; unfortunately, it requires an enormous amount of 
computing resources. 

The variational method has also been applied to a rib guide to determine 
the rib field profile. The accuracy of this technique is dependent upon the choice 
of the trial field function. Taking the trial function and varying its constants imtil 
an integral relation with the propagation constant is maximized, produces the field 
result. This is an application of the Rayleigh-Ritz procedure based on an expansion 
of a set of independent basis functions [30]. An integral relation is derived from the 
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scalar wave equation and defined as 

2 |$(a:,2/)|^da - J |Vt$(x,y)|^da 

Jmx,y)\'^da 

where the integration is taken over the guide cross section, $(rr,j/) is the field trial 
function, and is defined as ^ where A<, is the wavelength in free space. The 
variational method has received a fair eimoimt of attention recently [32], [33], [34]. 

The beam propagation method (BPM) is used to solve the scalar equation 
and cannot be adapted to solve the vector wave equation. Yevick and Hermansson 
compared the results of the finite element method with the beam propagation method 
[30]. Accurate results can be obtained for straight guides, but the BPM becomes 
inaccurate when applied to curved sections. The BPM is particularly useful for 
determining the coupling between two adjacent, parallel waveguides. 

The Fourier method has also been applied to a rib waveguide [35], [36], [37]. 
The electric field is calculated using a series of sine functions. Inserting this field 
expansion into the scalar wave equation and exploiting the orthogonality produces a 
matrix equation. Using sine basis functions, the boundary condition of zero field at 
infinity is not satisfied. Thus, the field must be enclosed in an artificial rectangular 
boundary, such that the accuracy of the results depends on its size [16] and the 
number of expansion terms used. 

After exploring the different methods for determining the field profile of a 
rib waveguide, three methods were chosen. The EIM and Marcatili’s approximation 
method were chosen because both have demonstrated accurate results for a rib waveg- 
uide. The Fourier method was chosen as the numerical technique because it produces 
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smaller matrix equations than the other numerical techniques [16]. This reduces the 
computing time required for the calculations, while maintaining accuracy. Due to 
the flexibility of the three chosen techniques, they can be applied to both coupling 
components. 

2.2.2 Fiber Field C^dculation 

The exact solution for the propagation characteristics of an elliptical-core 
fiber involves complicated transcendental dispersion equations. These equations re- 
quire numerical analysis to determine the fiber propagation constants. Consequently, 
most of the recent work on elliptical-core fibers has been based on approximation or 
numerical techniques [38]. 

Yeh [2] provided an exact analysis of an elliptical-fiber using Mathieu and 
modified Mathieu functions. From this analysis, a numerical technique is needed to 
produce the solutions. Using Yeh’s analysis, Dyott et. al. [39] developed a computer 
program to solve for the propagation characteristics of the fiber using the finite 
element method. This solution is very complex and is not conducive to the coupling 
efficiency calculation. 

The EIM is applicable to elliptical-core fibers [18], [40]. The solution is a set 
of spatially dependent analytic functions which are conducive to evaluation of the 
overlap integral. Fiber birefringence can also be calculated using this method. 

Another approximation method which can be applied to the fiber is Mar- 
catili’s technique [41], [42]. This method has been proven to accurately predict the 
propagation behavior of a highly elliptical-core fiber [43]. Birefringence results using 
this method compared favorably with those obtained by Dyott [39]. 
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The Fourier method has also been used to calculate the field profile of an 
elliptical-core fiber. This method utilizes an orthogonal set of basis functions to 
convert a differential equation into a set of simultaneous linear equations [44], [35], 
Results from this method accurately predicted the modal propagation constant and 
the field profile when compared with known results. 

From the above solution methods, the EIM, Marcatili’s approximation, and 
the Fourier method were chosen for this study. Both the EIM and Marcatili’s ap- 
proximation methods have been successfully applied to elliptical-core fibers. The 
numerical technique chosen was the Fourier method, based on its accurate results for 
elliptical fibers. Also, knowing a priori that (1) the Fourier method was chosen for 
the rib field calculation and (2) the overlap integral calculation was to be performed, 
choosing the Fourier method as the numerical technique for the fiber was a logical 
choice. All three methods produce field solutions conducive to calculating the overlap 
integral. 

2.2.3 Coupling Efficiency 

The goal of this research was to determine the feasibility of coupling light 
from an optical fiber into a rib waveguide. This required a comprehensive analysis of 
the propagation characteristics of both the elliptical-core fiber and rib waveguide over 
temperature. Results of this analysis were used to determine the theoretical coupling 
efficiency between the single-mode components for various alignment positions at two 
temperatures to determine interconnect feasibility. 

There are several methods of coupling a fiber to a waveguide, which include: 
prisms, gratings, evanescent field coupling, lenses, and end-fire coupling. All of these 
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methods were explored, but the most practical approach for an aircraft was end-fire 
coupling [45]. Prisms or gratings require the fiber to be incident at a particiJar angle, 
which is impossible in the vibrating environment of an aircraft enpne. Evanescent 
coupling requires the fiber to be stripped of its cladding and placed directly on top 
of the guide. Scratches on the upper surface of the guide cause scattering loss and 
precise alignment with the 2.7 fim wide rib is difficult. Fabricating a laser diode 
source directly on the device substrate [46] was not a viable option for this device 
since it was to operate in a temperature range from 20°C to 300°C, because the 
output wavelength of the laser diode would shift with temperature, corrupting the 
device output. End-fire coupling requires the fiber tip to be in physical contact with 
the guide cross section and precisely aligned, which is the most feasible coupling 
method for this application. 

Burns and Hocker [47] explored end-fire coupling in 1977, and discussed the 
issue of mode mismatch between components, which decreases the coupling efficiency. 
The coupling efficiency depends on the spatial overlap of the transverse fields of the 
optical modes in the fiber and channel guide, therefore any differences between their 
field profiles will produce a loss. The calculations performed by Bums and Hocker 
[47] approximated the fields of both components as Gaussian, to simplify the integral 
calculation. 

Approximating the fields of the two coupling components greatly simplifies 
the evaluation of the overlap integral. A closed form solution has been derived using 
the Gaussian approximation [48]. From this solution, the coupling efficiency between 
a GaAlAs laser diode and a LiNbOa device was determined. 
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Ghosh [49] performed an interesting theoretical probability analysis of the 
alignment between two components. Using a circular-core fiber which was smaller 
than the circular device, the probability of alignment was determined. For this case, 
100% coupling efficiency was possible because the fiber was smaller than the device. 
As long as the fiber was positioned entirely within the device cross-section, 100% 
coupling efficiency is achieved. This geometry was in contrast with the ceise tmder 
study where the fiber is larger than the device, being nearly two times as high. Thus, 
100% coupling efficiency can never be obtained for the device being studied. 

The majority of coupling research used LiNbOs for the waveguide, however, 
this material cannot be operated at high temperatures. Also, circular-core single 
mode fibers were used as the light input, thus simplifying the theoretical coupling 
calculation. Most papers simplified the coupling efficiency calculation by assuming 
the field of one or both of the components as Gaussian [50], [51]. 

Chung [52] addressed the problem of experimentally aligning the fiber with 
the guide, as well as the need for a polished waveguide end face to reduce scattering. 
The use of grooves fabricated directly on the waveguide for fiber placement improved 
alignment. A v-groove configuration is illustrated in Figure 2.1, where the angle 6 is 
the crystallographically defined angle of silicon. This approach is widely accepted and 
has been extensively used [53], [54]. Using v-grooves for fiber alignment is considered 
a self-alignment technique, because it reduces the necessary active alignment to one 
direction, z. Using LiNbOs devices, Murphy and Rice [55] experimentally fabricated 
a multiple fiber-waveguide alignment technique using a series of v-grooves on silicon. 
This design is suitable for both single-mode and multi-mode fibers. The required 
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groove width can be calculated using 

Wo = VSr — 2 A, (2.2) 

where Wo is the groove width, r is fiber radius and A is the distance from the silicon 
surface to the center of the fiber core. In practice, groove widths can be controlled 
to il fXTTi] this variation in the groove width produces a transverse offset causing 
coupling loss. The coupling efficiency was theoretically calculated assuming Gaussian 
mode profiles for both the fiber and waveguide. The fiber/guide system was thermally 
cycled between 20°C and 70°C, with minimal coupling efficiency change. 



Figure 2.1: V-groove in silicon for fiber alignment. 
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For alignment improvement, a fiber with a tapered diameter can be inserted 
in the groove and translated until the maximum coupling efficiency is achieved [56]. 
This scheme is illustrated in Figure 2.2. Using a LiNbOa device, coupling efficiencies 
of 70% or more have been experimentally achieved with this configuration [56]. 

Flip-chip alignment is similar to the v-groove concept. V-grooves in the 
substrate chip are used to align the fiber while another v-groove chip is placed on top 
and upside down such that the v-grooves are aligned. This method keeps the fiber(s) 
more securely in place. 



Figure 2.2: Silicon v-groove coupling structure 
with tapered fibers for improved alignment. 

Fabricating a lens on a fiber tip to improve field matching between the com- 
ponents is theoretically feasible, but experimentally difficult to achieve. Cai et. al. 
proposed the utilization of a coupling waveguide to improve field matching between 
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a single-mode fiber and a thin-film waveguide [57]. A 10 fj,m core polarization- 
maintaining fiber and 1.2 /im thick guide were used for coupling, which are much 
larger than the components being studied, with a calculated coupling efficiency of 
64%. 

More recent research involves gallium arsenide (GaAs) waveguide devices 
[54], [58]. Vezzetti and Munowitz [59] explored coupling between a single-mode circular- 
core fiber and a strip-loaded waveguide composed of GaAs using the Fourier method. 
They maximized the predicted coupling efficiency by altering the strip-loaded waveg- 
uide physical dimensions and indices so that its field more closely matched that of 
a circular-core fiber. Their research differed from this dissertation in that (1) the 
waveguide was a strip-loaded channel, (2) the fiber was circular-core, (3) the mate- 
rial used was GaAs, and (4) parameter changes with temperature were not explored. 

Experimentally, using a single-mode fiber, a 0.5 dB coupling loss to a GaAs 
device was reported, which corresponds to an 89% coupling efficiency [54]. This work 
was purely experimental, and no calculations of the fields of the components were 
performed. Another experimental group used a polarization-maintaining fiber and 
determined the coupling efficiency between it and a LiNbOa device. Intrinsic loss 
due to the field mismatch was approximately 0.5 dB [53]. 

A relevant paper compared the calculated coupling efficiency between four 
rib waveguides and a circular-core single-mode fiber using three different methods 
for the rib: the EIM, finite difference, and function fitting [60]. The function fitting 
method utilized the Rayleigh-Ritz procedure. The waveguide modes were assumed 
TE and propagating in the z direction, which are the same assumptions used in 
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this work. However, the fiber used in the paper was a circular-core single mode 
fiber having a spot size of 5 fim, with the field approximated as Gaussian. Results 
of the work indicated that the EIM was most accurate for ribs with small height 
and large widths, similar to the rib studied in this dissertation. Coupling efficiency 
variances due to transverse offsets and temperature change were not explored by 
these researchers. 

Only a few researchers addressed coupling a polarization-maintaining fiber to 
a waveguide. Of those, most approximated the field as Gaussian, and no calculations 
were performed to determine the coupling efficiency versus transverse offsets [61]. 

Thus, of the literature surveyed, none calculated the coupling efficiency using 
rigorous field solution techniques for both the fiber and rib waveguide to calculate 
overlap integral, and then determined changes in the coupling efficiency with trans- 
verse offsets at two temperatures. 

2.3 Ring Resonator 

The rib channel investigated for this research was the input to a micro- 
machined ring resonator pressure sensor, shown in Figure 2.3. There are several 
modulation techniques which can be used to deduce a measurand. These have been 
described elsewhere [62]. For the sensor investigated, wavelength modulation was 
used. Light propagates from the input rib channel waveguide through the first Y- 
branch and is split at the second Y^-branch. A portion of this light then travels 
around the ring where it crosses the diaphragm. The light continues around the ring 
to the first Y-branch and interferes with the incoming light. The effective index of 
the ring, rie, is altered when the diaphragm is perturbed by an applied pressure. The 
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round trip phase shift of the ring is given by [63] 

27T 

( 2 - 3 ) 

where A<, is the wavelength of the source and 71^ is integrated around the perimeter 
of the ring having length L. At resonance, the ring and input light constructively 
interfere such that = 2tti71j where m is an integer, which causes the transmissivity 
of the ring resonator to be a maximum. When (f) — 1 tx{th 4- Aj, the transmissivity is 
a minimum. This sensor exhibits cross-sensitivity to temperature which may be re- 
duced by shortening the path length [64] . Because the sensor is polarization sensitive, 
the input fiber must be polarization-maintaining. The ring resonator configuration 
can be used to measure other parameters by replacing the diaphragm in Figure 2.3 
with an appropriate sensing mechanism. 



Figure 2.3: Integrated optic ring resonator interferometer. 




CHAPTER III 


COUPLING COMPONENTS - RIB GUIDE AND ELLIPTICAL FIBER 

3.1 Introduction 

The channel waveguide discussed was the input stage of an integrated-optic 
pressure sensor. The pressure sensor was based on a strain sensitive ring resonator, 
combined with a micromachined silicon diaphragm. This sensor was polarization 
sensitive, thus interrogation with TE or TM polarized light each produced a different 
output [63]. In this investigation, only the incident TE mode was analyzed. The 
rib guide and elliptical-core optical fiber were the two components studied such that 
the theoretical coupling efficiency of the light output from the fiber into the input 
rib guide was calculated, with the goal being to maximize the coupling efficiency 
between the two. Both coupling components were single mode. The sensor design 
and fabrication was performed at the University of Cincinnati [63]. The physical 
parameters of both the sensor input and elliptical-core fiber, illustrated in Figure 
3.4, are discussed in the following sections. 

3.2 Rib Waveguide 
3.2.1 Materials 

Materials for the rib waveguide were judiciously chosen, consisting of a semi- 
conductor substrate followed by an isolation layer, guiding layer, and an overlayer. 
Each film layer is discussed in this order which was also the deposition order. 
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Figure 3.4: Cross sections of rib waveguide and elliptical-core optical fiber. 


3.2.2 Substrate and Isolation Layer 

Silicon was chosen as the substrate material due to its excellent waveguiding 
and mechanical properties [65]. For field confinement via total internal reflection, 
the index of the guiding region must be larger than the surrounding layers. The 
index of refraction of silicon at 830 nm is 3.67 [66] at room temperature. Because 
the index of all materials varies with wavelength, it is referenced to a wavelength . 
Unless otherwise indicated, all indices discussed in this dissertation are referenced 
to 830 nm. The index of silicon is higher than most candidate films used in high 
temperature applications. Therefore, an isolation layer of low index material must 
be deposited between the guiding layer and the substrate in order to utilize a silicon 
substrate and silicon micromachining technology. For this case, silicon oxynitride 
(SiON) was chosen as the guiding film, making silicon dioxide (Si02) a natural choice 
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for the isolation layer. Si 02 is dense, smooth, highly adherent to silicon substrates, 
uniform in refractive index and thickness, simple to produce, has a lower index than 
SiON, and can be fabricated nearly defect free. 

Because the isolation layer was finite, there was loss due to substrate cou- 
pling. Consequently, the waveguide’s modes were not strictly bound modes, but 
rather leaky modes [9], [67]. To minimize this effect, the thickness of the Si02 layer 
was chosen to nearly isolate the guide from the substrate. A theoretical analysis was 
performed [68] to determine the required isolation layer thickness to produce a 0.1 
dB/cm loss. According to the graphical results shown in Figure 3.5, the isolation 
layer thickness must be at least to 2.0 fim for a 0.516 fim thick guide to reduce the 
substrate coupling loss to 0.1 dB/cm. For the sensor, the isolation layer was 3.1 fim] 
thus, the guiding layer did not ‘see’ the substrate and the effects of substrate coupling 
were negligible. For modeling purposes, the isolation layer was considered to be the 
substrate. 

3.2.3 Guiding Layer 

Silicon nitride, SisN 4 , was first explored as the waveguiding core material 
because of its ability to withstand high temperatures without degrading. However, it 
was discovered that Si 3 N 4 films tend to crack for thicknesses exceeding approximately 
0.2 fim. Also, there is a large difference in the index of refraction between Si 3 N 4 
(n « 2.0) and an optical fiber (n « 1.5) which is not conducive for coupling, because 
Fresnel refiection would add a significant contribution to the loss. Using Si 3 N 4 , the 
calculated coupling efficiency would generally be less than 10%, which is unacceptable 
for this sensor application [69]. 
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Figure 3.5: SiOa isolation layer thickness vs. SiON guide thickness. 

To maximize the fiber-guide coupling efficiency, the cross sectional area of 
the two components should be nearly identical, as well as the core indices. Therefore, 
SiON [70], [71] was the material of choice for the guiding layer. The benefits of this 
material include the following: (1) film index can be varied from 1.46 - 2.02 with 
the fabrication process, (2) high quality deposited films, (3) ability to deposit thicker 
films as compared to Si 3 N 4 , and (4) better index matching to optical fibers. 

Ideally, the film thickness should equal that of the input fiber (1.3 (Jliti), 
but Low Pressure Chemical Vapor Deposited (LPCVD) SiON films deposited thicker 
than approximately 0.5 fim become multimode, and cannot be used for a ring res- 
onator sensor. The actual film thickness of the SiON layer was 0.516 /xm. Using 
photolithography, portions of the film were selectively etched to 0.426 fim to pro- 
duce the rib. This etch depth was chosen to provide optimal field confinement while 
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maintaining a single-mode guide. The rib width was chosen to optimize the sensor 
performance, taking into accoimt bends in the ring, and attempting to closely match 
the 2.6 fim major axis diameter of the fiber. For the device to remain single mode it 
was required that the channel width be approximately 2.4 fim. A mask was designed 
for this width, however, due to variations in the fabrication process, the actual width 
was 2.7 fxm. The thickness of the guide was chosen to match that of the elliptical 
fiber. In this study, the fiber diameters along the major and minor axes were 2.6 fJim 
and 1.3 fim, respectively. 

3.2.4 Overlayer 

An overlayer of LPC VD Si02 was deposited on top of the rib guide structure. 
This overlayer reduced the guide propagation loss by insulating the guiding film from 
the exposed surface, which minimized scattering effects due to surface imperfections. 
Similar to the isolation layer thickness, the overlayer thickness was chosen such that 
it appeared infinite to the guide. The effective index of a four layer guide consisting of 
(Si 02 \Si 0 N\Si 02 \air) was calculated for a Si 02 overlayer thickness varying between 
0 and 2.0 /xm. As illustrated in Figure 3.6, there was minimal change in the effective 
index for an overlayer oxide thickness exceeding approximately 0.5 /zm. Because the 
upper Si02 layer was 0.509 fj,m thick, it was modeled as an infinite layer. Thus, the 
rib guide structure was modeled as a three-layer structure as shown in Figure 3.7. 

3.2.5 Deposition Procedure 

The fabrication process, performed at the University of Cincinnati [63], be- 
gan with a 2” silicon wafer which was polished on both sides. The diaphragm was 
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Effective Index v. Overlayer Thickness 



Figure 3.6: Guide effective index calculated 
for various Si 02 overlayer thicknesses. 


fabricated first, then the wafer was stripped of any native oxide and placed in an oxi- 
dation furnace where a thermal oxide was grown in a wet oxygen environment. After 
the oxide was grown to the desired thickness, the wafer was placed in an LPCVD 
chamber to deposit the SiON guiding layer. 

Compared to other techniques, such as sputtering or Plasma Enhanced 
Chemical Vapor Deposition (PECVD), LPCVD produced the highest quality guiding 
films with the lowest loss [70], [71]. Increasing the amount of oxygen gas flow dur- 
ing deposition decreases the film refractive index, with the lower index limit being 
that of pure Si02, which is 1.453 [66]. Actual film refractive index tends to vary 
between runs even on the same equipment using identical conditions; hence there is 
no standard procedure to fabricate a film with a given thickness and refractive index. 




25 


3 Layer Rib Waveguide for Analysis 


Cladding Layer n = 1 .453 

Guiding Layer n = 1 .553 

Cladding Loyer n == 1 .453 

_ I I I I I I ^ 

- 2-10 1 2 

X (mic rons) 

Figure 3.7: Three layer structure for analysis. 

During the deposition run, the wafers were rotated to ensure uniformity of the film 
thickness and index. 

After the guide film was deposited, the ring resonator was fabricated using 
photolithography. Using a mask to define the rib sections, photoresist was deposited, 
and the rmprotected SiON regions were etched in acid by 0.9 /xm. After the rib was 
formed, an Si02 overlayer was deposited. Because Si02 films grow extremely slowly 
over SiON in an oxidation furnace, the Si02 overlayer was deposited using LPCVD. 

3.2.6 Parameter Variations with Temperature 

The effective index and field profile of the rib are expected to change with 
temperature. Three mechanisms can produce this change: (1) thermal expansion of 
the layers causing the width and thickness of the structure to vary, (2) stress induced 
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by the mismatch in the coefficient of thermal expansion of each layer, altering the 
refractive index via the photoelastic effect, and (3) index of refraction change with 
temperature. All of these mechanisms were taken into account to determine how the 
rib guide parameters changed with temperature [63]. 

For each layer of the rib waveguide, values for the thermal change in the 
refractive index, — » coefficient of thermal expansion Ethermai ’ 

photoelastic coefficients (pn and P 12 ), and the Poisson ratio (t/) were determined. 
These values can be foimd in the literature for all the rib constituents, except SiON. 
Values for silicon and Si02 vary with temperature, with averages over temperature 
given in Table 3.1 [72]. 

Compared to Si02 and crystalline silicon whose thermal properties are known, 
SiON is not strictly defined since its index can vary from 1.453 to 2.008 [66]. Its ther- 
mal properties vary as well and are dependent upon both the substrate material and 
film thickness [70]. The index of the SiON layer is 1.553 and its temperature sensi- 
tivity was approximated as that of Si 02 , whose refractive index is 1.453. 

For 0.1 fxm thick LPCVD SiON films on quartz (Si02), an experimental 
data set indicated that the film refractive index varied with temperature as 1.5 x 
10~‘^/°C, whereas for SiON films on a silicon substrate, it varied as —1.1 x 
[70]. This rib structure had a silicon substrate with a thick Si02 isolation layer, hence 
the refractive index dependence with temperature was expected to be between these 
two results, which that of the Si 02 satisfies. 

The following development is from De Brabander [63]. Since the rib layers 
were much thinner than the silicon substrate, changes in the film dimensions parallel 
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to the substrate {x and z directions) were the same as that of a substrate without 
any films, which is, 

^Xjthermal “ ^z^thermal “ 2.33 X 10 (^•^) 

If the thin Si02 and SiON films were removed from the substrate, they would expe- 
rience a strain which is expressed as 


^x, thermal — ^z, thermal — ^y,thermal — 0.55 X 10 


-6 


m 


m- °C' 


(3.5) 


Hence, there is a strain parallel to the substrate face which changes the refractive 
index via the photoelastic effect, denoted as 


^Xyphotoelastic — photoelastic — (2.33 0.55) X 10 


-6 


m 




(3.6) 


The substrate stretches or compresses the film and induces an additional strain in 
the y direction, through Poisson coupling, which is expressed as 


£y, photoelastic — ^ ^{,^x,photoelastic ^z, photoelastic') 


m- °U 


where u is the Poisson ratio, with value 0.17 [72]. 

For the rib waveguide, the thickness and index changes with temperature 
were small, hence applying each effect separately and adding produces the same 
result as applying all the effects simultaneously [63]. For example, the thickness 
change in the y direction is determined by the sum of the thermal and photoelastic 
strains, as given by 


^y, total ^y,thermal ^y, photoelastic — 1.791 X 10 


-7 


m 


m—°C 


(3.8) 
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For each parameter of interest, the expression and result for a \°C change of temper- 
ature is given in Table 3.2, where the abbreviation £f,pe where i = x,y,or z was used 
to denote Sx^'photoelastic’ 


3.3 Fiber 

3.3.1 Characteristics 

Compared to multimode fibers, single-mode fibers have much smaller core 
diameters. Thus, alignment is more critical when coupling a single-mode fiber to 
an integrated-optic device. Circular-core single-mode fibers can be characterized as 
having V < 2.40, where V is defined as 

V = (3.9) 

and a is the core radius [74]. Due to the symmetric geometry of a circular fiber, 
the field profile can be approximated as Gaussian. A circular-core single-mode fiber 
actually propagates two orthogonal fundamental modes, corresponding to TE and 
TM polarization [39]. Thus, to maintain the polarization state of the input light 
through a fiber requires an alteration of the standard circular fiber. Preliminary 
calculations were conducted using the specification data from commercially available 
circular-core single-mode fibers. For a circular-core radius of 4 /xm, the calculated 
coupling efficiency was unacceptably low, 45% [75]. In an effort to improve the 
coupling efficiency, other fibers were explored. 
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An elliptical-core polarization-maintaining fiber was chosen because it was 
expected to provide the highest coupling efficiency due to its geometry and refractive 
index being closely matched to the guide’s. Polarization-maintaining fibers can be 
one of two types: (1) a circular-core fiber under stress from the cladding or (2) an 
elliptical-core fiber. For this study, an elliptical-core fiber was used because its di- 
mensions and field profile more closely matched that of the guide, reducing coupling 
loss. These types of fibers can be used in coherent optical communication systems, 
fiber-optic sensors, and for coupling into polarization-sensitive integrated optic de- 
vices [27]. The latter case applies to this research. Elliptical-core fibers have received 
more attention recently due to the availability of longer wavelength sources and de- 
tectors. Prior to this, the lower wavelengths available required the elliptical-core 
dimensions for polarization preservation to be too small for fabrication [39]. 

An elliptical core-fiber has two preferred directions of polarization; along 
the major axis and along the minor axis [74]. Polarized light coupled into the fiber 
will propagate as two orthogonal modes, denoted EH and HE. EH denotes modes 
with strong axial magnetic field strength, H^, corresponding to a quasi-TE mode 
whereas HE denotes modes with strong axial electric field strength, E^, correspond- 
ing to a quasi-TM mode [76]. Each mode travels at a different velocity, known as 
birefringence. The fiber modal birefringence is defined as [74] 

n ~ /^y) 

^ 2 ^ ■ 

where j3x is the TE mode propagation constant and /3y is the TM. Because Px > Py, 
the EH mode is more strongly guided [39]. 
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By deliberately inducing a large degree of birefringence in a fiber, establishes 
a preferred polarization direction into which all the power is launched [77]. The beat 
length of a fiber is inversely related to its birefringence, and is defined by 



For short beat lengths or high birefringence, transfer of power to the orthogonal mode 
is precluded. Hence, the output remains polarized in the direction of the preferred 
axis. 

Highly birefringent fibers are less susceptible to mode transitions, such that 
the output polarization mode is the same as the input polarization. Less birefrin- 
gent fibers are susceptible to perturbations in fiber dimensions and indices along its 
length such that an input TE polarization may transfer to the TM polarization state, 
which is an undesirable property for this application. Highly birefringent fibers are 
considered polarization-maintaining [74], the desired property for this application. 

3.3.2 Physical Parameters 

The fiber used was commercially available [78]. Data from the manufacturer 
indicated that the diameter of the major and minor axes of the core ellipse were 2.6 
fim and 1.3 /zm, respectively. The nominal values of the core and cladding indices 
were 1.484 and 1.450, respectively. The numerical aperture (NA) of the fiber, defined 
as 

NA = ^nf^, - 


(3.12) 
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was 0.316. These parameters will vary slightly between batches, but for this study, 
the nominal values were used. A fiber is characterized as “weakly guiding” if the 
difference between the core and cladding index is small. For this fiber, the difference 
was approximately 2.3%, hence it can be classified as weakly guiding [76]. 

3.3.3 Parameter Variations with Temperature 

The elliptical-core fiber investigated had a silica cladding with a germanium 
doped core. Generally, the cladding material has a slightly higher coefficient of 
thermal expansion than the core [79]. However, lacking data for the doped core, 
the same coefficient of thermal expansion for the cladding was used, 5.9 xl0~^— 
which is valid from 0°C to 300°C [78]. Thus, the stress induced by differences in 
expansion between the core and cladding was neglected. Similarly, the core change 
in refractive index due to temperature was approximated as that of the cladding. 
Parameter changes for a 1°C change in temperature are given in Table 3.3. 

3.4 Summary 

Both the rib waveguide and polarization-maintaining fiber were discussed in 
this chapter. Reasons for choosing the material for each rib layer and its respective 
dimensions were given. The elliptical-core fiber was selected due to its geometry 
closely matching the rib and its polarization-maintaining properties. The index and 
dimensional variations of both components with respect to temperature were dis- 
cussed. Results of these calculations will be presented in Chapter V. The methods 
used to independently calculate the electric field profile and propagation constants 
of both coupling components will be explored in the next chapter. 
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Table 3.1: Temperature Dependence of Silicon and Silicon Dioxide 


MATERIAL 


^thermal [65] 

Pii [72] 

Pi2 [72] 

v[n] 

Silicon 

1.5 X lO-'^ 

2.33 

NA 

NA 

NA 

Si02 

1.28 X 10-^ 

0.55 

0.121 

0.27 

0.17 


Table 3.2: Parameter Change Calculation 


PARAMETER 

Expression For Temperature Dependence 

EBH 

W 

^x, total 


toverlayer 

^y, total 


iunetched 

^y^total 

-9.241 X 10"“ iim 

^etched 

^y^total 

-7.632 X 10“* iim 

'^overlayer 

“(^^eWayer/ 2 ) {^x.pePn + [^z,pe + ^y,pe]Pl 2 } + 5 ^ 

1.2035 X 10-® 

T^rib 

~{'^rib/^) {^x.pelPu ^z,pe +^y,pe] P12} dv 


'^isolation 

~{'^%olatian/*^) {^x,pePu i^^,pe + ^y,pe] P12} + 



Table 3.3: Parameter Change Calculation for Fiber 


PARAMETER 

Temperature Dependence 

Calculated change 

a 

^Xy fiber 

7.67 X 10~“ fj,m 

b 

^y, fiber 

3.84 X 10“^ fi,m 

'^core 

dn 

Jt_ 

1.28 X 10-^ 

'^cladding 

dn 

Jr 

1.28 X 10-^ 



































CHAPTER IV 


METHODS USED TO DETERMINE THE ELECTRIC FIELD 

PROFILES 


4.1 Introduction 

As fiber optic and integrated-optic device technology areas have developed, 
so has the need for understanding the modal characteristics of these structures. This 
need arises from the desire to accurately predict the performance of a device prior 
to fabrication. In particular, rib, channel, and strip-loaded waveguide devices have 
received more attention due to improvements in fabrication, particularly microma- 
chining. These structures have been developed for practical uses as sensors, couplers, 
and electro-optic devices. 

4.1.1 Three-Layer Planar Waveguide - Analytical Solution 

To provide a basis for understanding a rib waveguide, a discussion of the 
three-layer planar guide solution will be given following the development of Marcuse 
[10]. Each dielectric layer is characterized by its index of refraction, n. The index 
of refraction is a dimensionless quantity relating the speed of light, c, to the phase 
velocity, v, of the electromagnetic fields within the dielectric material. That is, 

n = ~ = \/ fJir^ry (^•^3) 

where e = EoEt and /i = /Xo/^r are the permittivity and permeability of the medium 
with So and fio denoting the free-space values. The three-layer planar waveguide 
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is shown in Figure 4.8 where n\ is the core index of refraction and 7i2 and nz are 
the cladding indices. The x, y, and z coordinates illustrated in the figure are used 
consistently throughout this dissertation. 

To confine the electromagnetic field within the guiding layer requires that 
Til > Ti2 and Til > Generally, the overlayer has a smaller index than the lower 
cladding such that the confining condition can be written as tii> Ti2> nz. A planar 
waveguide confines the field only in the y direction while the field in the x direction 
extends to infinity. The field propagates in the z direction via total-internal-reflection 
at each of the dielectric interfaces, ni/ ri2 and ni/ nz interfaces. 



Figure 4 . 8 : Three-layer planar waveguide. 
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The derivation begins with Maxwell’s equations for linear, homogeneous, 
and isotropic nonmagnetic media, 

5E 

V X H = (4.14) 

da 

VxE = (4.15) 

where H and E are the time-dependent magnetic and electric field vectors. The anal- 
ysis is simplified by assuming no field variation in the x direction, denoted d/dx = 0. 
For brevity, only the transverse electric (TE) development for guided mode solutions 
will be given here. The TM development is analogous to the TE, hence only the TM 
results will be given. 


4.1.2 TE Mode Development — Guided Modes 

TE modes have only three field components, Ej., Hj,, and H^. When E„ = 
Ez = Hj; = 0, equations (4.14) and (4.15) become 


dHy dH, 

dz dy ~ dt ’ 

dE^ dHy 

dz dt ’ 

dE^ _ dH, 

dy ~ dt' 


(4.16) 

(4.17) 

(4.18) 


The time and z dependence of E and H are given by which describes a wave 

traveling in the positive z direction with phase velocity u = uj/P, where lo = 2-nf is 
the radian frequency of the wave and P is the propagation constant. When 
is inserted into equations (4.16) - (4.18), they become 


dHz . 2n 

-tPHy - = IUEoU^Ex, 


(4.19) 
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i0Ex 

dEx 

dy 


-iUfloHy, 

-iwfioHz. 


(4.20) 

(4.21) 


Rewriting equations (4.19) - (4.21) in terms of Ex yields 


p ^ S_„ I i SH, 

* uSoTi^ ^ u€on^ dy ’ 


Hy = 


Hx 


Xi 


--^E. 

Uflo 
i dEx 


(4.22) 

(4.23) 

(4.24) 


ufio dy 

Substituting equations (4.23) and (4.24) into equation (4.22) and multiplying by iw/Zo 
produces 

-P^Ex + = -uPeolJion^Ex, (4.25) 

which upon rearrangement is the one-dimensional wave equation, 


^ + {n^kl-P^)Ex = ^. 

where = (27 t/Ao)^, A<, is the wavelength of light 

defined as 


(4.26) 


in free space, and n is 


n = riz 

for (y > 0), 

(4.27) 

n = ni 

for (0 > y > — d). 

(4.28) 

n = 722 

for (y < —d). 

(4.29) 


Solutions to equation (4.26) must satisfy the boundary conditions at the dielectric 
interfaces, y = 0 and y = —d. The tangential E and H fields must be continuous 
at the dielectric discontinuities. That is. Ex and Hz are continuous at y = 0 and 
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y — —d and the solutions for Ex must approach zero as y approaches ± oo. Thus, 
solutions to Ex with implicit time and z dependence, are of the form 


E.{y) = 

Ae~^v 

for y > 0, 

(4.30) 

E.{y) = 

[i4cos(/cy) + Bsin{Ky)] 

for 0 > y > —d, 

(4.31) 

EAv) = 

[/lcos(Kd) — Bsin{Kd)] 

for y < —d. 

(4.32) 


These are the solutions of Maxwell’s equations producing the electric field of an 
antisymmetric three-layer planar waveguide in the three dielectric regions. If nz does 
not equal the guide is antisymmetric; conversely, if nz equals TI2, the planar guide 
is symmetric. The parameters 5, k, and 7 are given by 


= yJP^-nlkl 

(4.33) 


(4.34) 


(4.35) 


To solve equation (4.26), the boundary condition for Hz must be enforced; hence Hz 
must be determined in the three dielectric regions. The results for Hz in regions 713, 
ni, and ri2, with implicit time and z dependence, are as follows; 


H.(y) 

H,{y) 

H.iy) 


[i4sin(K7/) — Bcos{Ky )] , 


U/Xo 

ia 


[i4cos(/cd) — Bsin(/«i)] 

UJflo 


(4.36) 

(4.37) 

(4.38) 
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4.1,3 TE Results for a Planar Waveguide - Guided Modes 

Imposing the continuity requirement of Hz at the two interfaces produces 
two equations 


6A + kB = 0, (4.39) 

[/csin(Acd) — cos(/cd)]i4+ [Kcos{Kd) + sin{Kd)]B = 0. (4.40) 

Solving this system of equations yields the eigenvalue equation for a three-layer planar 
waveguide, 

> ( 4 - 41 ) 

where 6, 7 , and k are functions of the propagation constant P, as given in equations 
(4.33) - (4.35). 

Solutions to this transcendental equation, equation (4.41), can be exactly 
determined. The number of solutions obtained yields the number of guided modes 
which can be excited in the guide. For a single-mode three-layer planar guide, there 
is only one solution, denoted the fundamental mode. For many applications it is 
desirable for the device to operate as single mode because higher order modes have 
larger attenuations. 

For the symmetric slab waveguide case, = U 2 and 7 = <5, and equation 
(4,41) becomes 



(4.42) 
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Equation (4.42) is second-order in tan(Kd/2), yielding two solutions, 

tan(/cd/2) = 7//C, (4.43) 

and 

tan(Kd/2) = -/c/7. (4.44) 

Equations (4.43) and (4.44) correspond to the even and odd modes of a planar waveg- 
uide, respectively. The fundamental mode is denoted mode zero, the next higher 
mode is one (odd), followed by two (even), etc. Generally, the thicker the guide 
region, the more modes it will support. These equations yield exact solutions for /3 
and the effective index of the guide is then calculated using 

He = P/ko- (4.45) 

4.1.4 TM Results - Guided Modes 

For the transverse magnetic (TM) incident case, the field has Hx, Ey, and Ez 

components. Using the same procedure as in section 4.1.1, the eigenvalue equation 

for the asymmetric planar guide is given by 

nj«(ri37 -f- 7126 ) 

{n^nzKf' — n\'^8 

For the symmetric slab case, ri2 =nz, the eigenvalue equation can be determined for 
the even and odd modes. Even modes are given by 

tan(/cd/2) = (4.47) 

and odd modes by 

tan(/cd/2) = - (4.48) 
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Boundary conditions require Hx and to be continuous at y = 0 and y = —d and 
Hx to vanish as y approaches ± oo. The electric field profile can then be determined 
from the Hx solution. 

4.1.5 Radiation Modes 

Radiation modes are solutions to Maxwell’s equations that also satisfy the 
boundary conditions. The radiation field is considered a mode not confined to the 
core but rather one that extends to infinity in the y direction. For this type of mode, 
the propagation constants are not restricted to a discrete set of values, because they 
are related to the plane wave angle of incidence. These propagation constants form 
a continuum, hence the radiation modes do so as well. 

The radiation modes can be further divided into substrate radiation modes, 
air radiation modes, and evanescent modes. All of the modes can be categorized by 
the value of the propagation constant. For guided modes, the propagation constant 
is in range kni < /? < ku 2 . For substrate radiation modes, it is confined to ku 2 < 
(5 < kns, for air radiation modes it is confined to 0 < j3 < kris, and evanescent modes 
are denoted by imaginary f5 values [81]. Because this research was concerned with 
the guided field profile, radiation modes were discussed for completeness but were 
irrelevant to the coupling efficiency calculation [10] . 

4.1.6 Rib Waveguide 

For this study, a rib waveguide structure was investigated. A rib waveguide, 
illustrated in Figure 4.9, is a planar-slab waveguide with a raised height section of 
given width [29], which provides lateral confinement. Thus, this structure confines 
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the light both vertically and horizontally, analogous to a channel guide, illustrated 
in Figure 4.10. That channel can be readily analyzed due to its simple geometry; 
however, rough sidewall edges are inevitable when fabricating such a device. These 
rough edges greatly increase the guide’s propagation loss, making it imdesirable for 
aircraft applications since the higher the propagation loss [80], the more input light 
power is required for interrogation. Hence, the rectangular channel guide works well 
in theory, but cannot be fabricated with a low enough propagation loss for practical 
use. 

Experimental results using various channel geometries such as channel, strip- 
loaded, rib, and buried rib, indicate that the buried rib waveguide produces the lowest 
propagation loss, which is desirable for a sensing device [63], [82]. A buried rib is a 
rib guide with an overlayer as illustrated in Figure 4.11. Propagation loss can be 
caused from leakage into the substrate, scattering at each guide/cladding interface 
and/or scattering due to imperfections in the film. Using the buried design decreases 
the field amplitude at the surface, thereby reducing scattering effects due to surface 
contamination, scratches, and dust. 

Unlike the planar guide, there is no exact solution to the two-dimensional 
Maxwell wave equation, due to the complex geometry of the rib channel waveguide. 
Therefore, approximation and numerical methods must be utilized. Several review 
papers have been written describing these methods [19], [20]. A discussion of the three 
chosen methods will follow in the next section. 
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Figure 4.9: Rib waveguide. 

4.1.7 Elliptical-Core Fiber 

The propagation constant of an elliptical-core fiber can be solved analytically 
using a series of Mathieu functions [2]. As such, the field components inside and 
outside the core must be represented by an infinite product of terms of Mathieu and 
modified Mathieu functions and the propagation constant determined from the roots 
of an infinite determinant [2]. A niunerical technique such as the method of successive 
approximations must be utilized to determine these roots. This numerical procedure 
is complex and time-consuming and the results are not conducive to the coupling 
efficiency calculation in this case. Thus, the same approximation and numerical 
techniques chosen to solve the propagation and field profile of the rib guide have also 
been used for the elliptical fiber. 
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4.2 Methods 

For the rib channel guide and optical fiber, three methods of calculating 
the propagation constant and electric field profiles were utilized: the effective-index 
method (EIM), Marcatili’s technique, and the Fourier method. All solve the scalar- 
wave equation. For completeness, a development of the wave equation, based on that 
of Marcuse [16], is given. Beginning with Maxwell’s equations, 


V X H = iwSon^E, 

(4.49) 

V X E = — iw/ioH, 

(4.50) 


the curl is taken of equation (4.50) and the result of equation (4.49) substituted on 
the right hand side, producing 


V X (V X E) = 


(4.51) 
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Figure 4.11: Buried rib waveguide structure. 


Using the relations 

V X (V X E) = V (V • E) - V^E (4.52) 

and 

V • (n^E) = E • Vn^ + n^V • E = 0, (4.53) 

one can obtain the vector wave equation for the electric field, which is 

V^E + V ^E . ^ j + n^fc^E = 0. (4.54) 

The guided mode solutions have a z dependence of e~*^; thus the Laplace operator 
is of the form 

The wave equation, equation (4.54), can be decomposed into a;, t/, and z components. 
Using the previously mentioned property that the field components can be expressed 
in terms of Ex and Ey, only two of the three equations are needed. Using this and 
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the assumption that 


dn - 


(4.56) 


produces the coupled equations for the x and y components of the vector wave equa- 
tion, which are 


dx"^ 

and 


+ 




+ [n^(rc, y)kl - ^0^] £?* + 2 


_5_ 

dx 


d^Ey d'^Ey r 2 / m 2 /d 21 n n ^ 


P ^ln(n) ain(n) 


p ^ln(n) 51n(n) 

r\ * ^y' 


= 0, (4.57) 


dx 


dy 


= 0. (4.58) 


dx“^ ■ dy^ ■ L ^ 

If the index difference between the core and cladding of the guide is small, 
such that 

A = « 1, (4.59) 

rii 

where ni is the core and 7i2 is the cladding index, then the component is considered to 
be weakly guiding [40]. For the rib waveguide A = 0.064 and for the fiber A = 0.023, 
which are both much less than unity, satisfying the weakly guiding condition. Thus, 
the last bracketed term in both equations (4.57) and (4.58) can be neglected. This 
yields the scalar- wave equation for Ex and Ey, 

+ [nHx,y)kl - 0^] EA..V) = 0, (4.60) 

and 


+ p(x, y)kl - E,(.,y) = 0. (4.61) 

The solution to these two scalar wave equations are identical. For quasi-TE modes 
propagating in the z direction, the electric field is of the form 


E(i,y,z) =i E,(x,y) e’*"‘ 


(4.62) 
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where Ex{x^y) is the field profile. The electric field of a TE mode has x and y 
components, however, the quasi-TE mode electric field has only an x component. 
All three of the chosen techniques, the EIM, Marcatili, and the Fourier method, 
produce solutions to the two-dimensional scalar wave equation (4.60), where n(x, y) 
is the spatially dependent function for the index of refraction and ko was previously 
defined. 

4.2.1 Effective-Index Method (EIM) 

Rib Waveguide 

The EIM reduces the two-dimensional scalar wave equation, equation (4.60), 
into two one-dimensional problems. This is achieved by first dividing the rib guide 
into three separate planar guides oriented perpendicular to the y axis, as shown 
in Figure 4.12. Slabs A and C, which have identical characteristics, correspond 
to the etched portion of the waveguide, while Slab B corresponds to the unetched 
portion. The effective index of each slab was determined using the procedure for 
the three-layer planar guide discussed in section 4.1.1 assuming TE incident light. 
Ni is the effective index of Slab B and Nji is the effective index of Slabs A and C, 
where Njjj = To determine the degree of confinement in the x direction, a 

symmetric slab waveguide, Slab D, is constructed perpendicular to the x axis. Due 
to the orientation of the guide, the effective index of Slab D is determined for the 
TM case. Slab D uses Ni and Nn as the core and cladding indices, respectively, and 
the rib width W as the core thickness. The analytic solution of Slab D yields the 
effective index, of the rib guide. 
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Figure 4.12: EIM constituent slab waveguides. 


While not applicable to this particular rib, it should be noted that when 
analyzing the slabs, it is possible for Slab B to be single mode while Slabs A and C 
are cutoff. If either Slab A or C is cutoff, there is no guided-mode solution for these 
planar regions. An approach to circumvent this situation is to assume an effective 
index of unity in these regions [83]. 

Quantitatively, for the TE mode, the transverse field is given by 

E* = ExxX + ExyV, (4.63) 


where Exx and Exy are the x and y components of the field profile. Substituting 
equation (4.63) into the vector wave equation, (4.54), produces a pair of coupled 
equations. 


+ inX - 0l)E^ 
+ {nX - 


r E„a(n^) 1 _ [ E^a(n^) ' 

dx dx dx ri^ dy 

d Exx d{n^) _ Exy 5(n^) 
dy dx dy dy 


(4.64) 

(4.65) 
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For Exx » Exy, by an order of A given in equation (4.59), the right hand side of 
these coupled vector equations is near-zero for small index discontinuities. These 
coupled vector equations can then be approximated by the scalar wave equation, 
equation (4.60). 

For a rectangular waveguide, the EIM can first be applied to either the 
rectangle oriented perpendicular to the x axis (short guide height), or the rectangle 
oriented perpendicular to the y axis (large guide height). Kumar et. al. [84] have 
compared the results for the same guide calculated both ways. They determined that 
the EIM produces more accurate results if applied to the shorter dimension first. 

The EIM uses the method of separation of variables to solve the scalar-wave 
equation, where the index profile is approximated as 

ri^{x,y) = n'^{x) + n"^{y). (4.66) 

As discussed for the rib, the 'field profile, Ex, is defined in equation (4.77). Using 
equation (4.66), the two-dimensional scalar wave equation can then be written in the 
form 

For the planar waveguide analyzed, n"{y) is denoted as 


n"{y) 

= ni 

for y < |6*j 

(4.68) 

n{y) 

= ri2 

for y < —b' 

(4.69) 

n"{y) 

= n3 

for y > b' 

(4.70) 
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and Y(y) satisfies the one-dimensional wave equation, 


1 (fY{y) 


+ 




= 0 . 


Y{y) dy^ 

Subtracting equation (4.71) from equation (4.67) produces 


(4.71) 


1 I /, 

X^x) dx^ \ 


n\x) + 1 


-y8M = 0, 


(4.72) 


where [3 is the propagation constant of the fundamental mode of the planar guide 
whose index profile is given by the bracketed term, 




PI PI 


K 

2 


kl 




for |a;| < a 
for Ixl > a 


(4.73) 

(4.74) 


which provides the x dependent index profile. 


n"^{x) = 0 


n'^{x) = — 


kl ^ 


for jrr| < a 
for |a:j > a 


(4.75) 

(4.76) 


Substituting the results for n"^{x) and n"^[y) into equation (4.66) produces the sep- 
arable refractive index profile, n3{x,y), the EIM uses to solve the two-dimensional 
scalar wave equation. This index profile is shown in Figure 4.13 and is larger than 
the actual index profile when |rr| > a and |?/| < b',hy n\- and smaller in the 

corner regions by — 714 . Since the field magnitude is smaller in the comer regions 
ko 

than the cladding regions, the effective index calculated using the EIM is always 
larger than the actual effective index, hence, this method consistently overestimates 
the propagation constant of the guide [84]. 
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Figure 4.13: EIM separable refractive index 
profile shown in its respective region - shorter 
dimension analyzed first. 


Following the same procedure, but starting with the rc-slab (larger dimen- 
sion) first, the EIM approximates the waveguide shown in Figure 4.14. In this case 
the index in regions y < —b' and y > b' is larger by nf - ^ and lower in the corner 

regions by ^ - n|. The error in the propagation constant is larger due to the overes- 
timation of the index in the cladding regions along the shorter rectangle dimension, 
where the field magnitude is greater than in the outer cladding region. Thus, when 
analyzing a rectangular guide oriented perpendicular to the x-axis using the EIM, 
the results are more accurate if the y-slab is analyzed first. 

The electric field profile of the rib is the product of the field distribution of 
Slab B, AT(x), and that of D, Y{y), i.e.. 


E,(x,y) = X{x) Y{y). 


(4.77) 
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The field function, X{x), is determined using equations (4.30) - (4.32) and Y{y) is 
determined from the eigenvalue equation (4.47), applying the boimdary conditions 
for the TM case. 

The EIM has intuitive appeal; however, because the guide is analyzed as 
three separate slabs, the field continuity conditions at the interfaces are not satisfied. 
Therefore, the results from this method are most accurate for ribs with small height 
and large width [60]. This condition was satisfied in this study as the guide etch 
depth was 0.09 fim which is small when compared to its 2.7 fim width. Several 
researchers have attempted to improve this method but none significantly change the 
results [22], [23], [24], [25], [26], [28], [85] and they tend to apply to specific cases. 


Fiber Field - EIM 

The fiber shown in Figure 4.15 is a single-mode polarization-maintaining 
elliptical-core fiber. Thus, the Gaussian approximation commonly used for circular- 


core fibers was not valid for this case [18], [40]. The electric field and propagation 
constant of the elliptical-core fiber were calculated using the EIM [18], where the fiber 
was approximated as a rectangle with dimensions a and b' such that the rectangle 
and ellipse have the same area and aspect ratio. Since an ellipse has an area of Tvab, 
these dimensions are given by 

/ y/n 

a = a-^ 

2 



and 



(4.79) 


Using the procedure given in section 4.2.2, the fiber field and effective index were 
determined for this rectangle. Marcatili’s method was designed to evaluate rectangles, 
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hence, for comparison purposes between the two approximation techniques, the fiber 
was modelled as a rectangle for the EIM. 



Figure 4.15: Elliptical-core fiber with pseudo-guide (rectangle). 


4.2.2 Marcatili’s Method 
Rib Waveguide 

For modes far from cutoff, the majority of the field energy is confined to the 
central core region. If this assumption applies, a pseudo-rectangular waveguide can 
be utilized to approximate the rib structure [15]. The cross section of the waveguide 
is divided into nine regions as shown in Figure 4.16, in which t is the etch depth and 
r is the thickness of the etched SiON layer, such that the unetched height is r t. 
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Figure 4.16: Marcatili’s method: pseudo- 

guide used to approximate rib guide. 


The pseudo-guide is assigned a refractive index profile which approximates 
that of the actual rib. This profile provides readier analysis because n^{x,y) is sepa- 
rable in the x and y coordinates, as defined by [15] 

y) = -I- n"^(y) - n\, (4.80) 

where 



n'(x) = 

ni 

IV 

for |a:| < y 

(4.81) 


n'{x) = 

U2 

W 

for |a:| > y 

(4,82) 

n"{y) 

= ni 

for 

-t<y <r, 

(4.83) 


= Tl2 

for (y < —t) or (y > r). 

(4.84) 


and 
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Analogous to the EIM, this method reduces the two-dimensional wave equa- 
tion into two separate one-dimensional problems. Substituting equations (4.81) - 
(4.84) into equation (4.80) produces an index, n{x,y), which is separable in x and y. 
Using this index distribution in equation (4.60) yields 

■^-^kln^{x)-/3l X{x)=0, (4.85) 

and 

^ + kln'^iy) - Pi y(y) = o, (4.86) 

where 

0i+Pl-kln\=p\ (4.87) 

When the mode is well-guided, the majority of the field energy is in the 
central core (ni region) which is the key assumption of this theory [10]. The refractive- 
index profile of the rib waveguide and pseudo guide differs in the outer shaded regions 
and sections D and H of Figure 4.16. The index of refraction in the shaded regions 
is — n\. In regions D and H the actual index is 

n = ri 2 for {-t < y < 0) (4.88) 

n = ni for (0 < t/ < r). (4.89) 

However, for this approximation, it is assumed that is the index throughout re- 
gions D and H. Because the modal power in these outer regions is small, the mode 
field profile of the pseudo-guide should closely resemble that of the rib guide. This 
assumption is generally applicable for ribs with large etch depths (t » r). How- 
ever, for this channel guide the etch depth (t) is small, but because the rib is wide 
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iyV » r + 1), it is expected that the pseudo-guide will provide an accurate field 
profile. 

Following the development of Varshney and Kumar [15], the field is given by 
equation (4.77), where the x dependent solutions are 

X{x) = Aicos(2fii-^ -9^ for|x|<y, (4.90) 

X{x) = ^ 2 «xp for |x| > y , (4.91) 

where 

(4.92) 

W y 

f^2 = (4.93) 

and 9 = 0 for a mode S 3 nnmetric (antisymmetric) in x. The y dependent 

solutions are given by 


y{y) 

- At, exp 1 

[=r) 

for y > r, 

(4.94) 

y(y) 

= Ac cos ^ 

:¥)+"-K¥) 

for — t < y < r. 

(4.95) 

y(v) 

= Aa exp 

(¥) 

for y < —t, 

(4.96) 


where 


76 = 

(4.97) 

7c = — ^2) 

(4.98) 

7a = 76. 

(4.99) 
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Assuming the dominant mode is quasi-TE, the boundary conditions are, 


n 


, , , dEJx,y) . W 

,Ex{x,y) and continuous at x = ± — , 

OX 2t 


(4.100) 


and 


Ex{x,y) and 


dExjx, y) 
dy 


continuous y = —t, r. 


(4.101) 


These boundary conditions are applied and an eigenvalue equation determined for 
the X and y dependent field functions, given by 


tan ^ 


nl Hi 


-//i + (p-l)| = 0, 


(4.102) 


and 


tan 


-1 


la 

7c 


+ tan 


-1 


76 

7c 


- 7c (l + 0 + (g - 1)7T = 0, 


(4.103) 


for the X and y dependent field functions, respectively. For the fundamental mode, 
p = q = l\n equations (4.102) and (4.103). Constants Ai, A 2 , A„, As, Ac, and Be 
are determined by satisfying the boundary conditions given in equations (4.100) and 
(4.101) and the electric field profile is calculated using the solutions for Pi and ^ 2 - 
Perturbation Correction - Guide Effective Index 
The error caused by approximating the rib guide with a pseudo guide can be 
reduced using a first-order perturbation correction to the propagation constant [14]. 
The first-order perturbation correction minimizes the difference between the actual 
guide and the pseudo guide. 
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The following development is from Varshney [15]. The perturbation cor- 
rection to the effective index, calculated using Marcatili’s method, is expressed as 


[15] 


= 


klJJ \Ex{x,y)^ S'n? dx dy 

j J \Ex{x,y)\'^ dx dy 


(4.104) 


where the integration is over all space and Sn^ is given by 

6n^ = {n\ —TI2) shaded regions (C), (E), (G), and (I), (4.105) 

6n^ = 0 otherwise. (4.106) 

Inserting P ' into the following equation produces the corrected value for /3, 

PI + P'^ = P\ (4.107) 

where Po is the propagation constant determined using Marcatili’s method. Because 
{n{ — rij) is small in this case, the correction term is small as well. Once P is deter- 
mined, the effective index of the guide is calculated using equation (4.45). 
Marcatili - Fiber Electric Field 

Using Marcatili’s technique to calculate the field profile of an elliptical-core 
fiber has been discussed in the literature with accurate results [14], [42], [43]. For this 
method, the elliptical core is approximated by a rectangle, for which an analytical 
solution can be found. The same rectangle dimensions chosen for the EIM are used 
for this method [42]. Using this rectangle as the pseudo-guide, Marcatili’s method 
as described in section 4.2.2, is used to determine the fiber field profile. The field is 
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given by equations (4.82) and (4.95) - (4.104), replacing t and r with b' , W with 2a', 
and Be — 0. 

Perturbation Correction - Fiber Effective Index 
To determine the perturbation correction to the propagation constant, equa- 
tion (4.104) is used. However, for the fiber, 6n^ is of the form [42], [43] 


6n^ = {n\ — Tij) 

for regions 1, 3, and 4, 

(4.108) 

6n^ = {n\ — n\) 

for region 2, 

(4.109) 

to 

II 

O 

otherwise. 

(4.110) 


where the regions are defined in Figure 4.17. Again, this calculation for the effective 
index is used only for comparison purposes between the methods. 



Figure 4.17: Elliptical-core fiber with re- 

gions defined for perturbation correction to 
the propagation constant. 
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4.2.3 The Fourier Method 

Rib Waveguide 

The scalar-wave solution applies to waveguides having small index discon- 
tinuities in the region of field confinement, which is the case being studied. If the 
waveguide had a large index difference between the core and cladding where the field 
intensity was large, a vector solution would have been necessary. Using Maxwell’s 
equations, elimination of either the magnetic or electric field vectors as shown in 
Chapter III can be achieved, to obtain the wave equations for either the electric 
field or magnetic field. Analogous to the procedure described in section 4.1.1, the 
transverse components of the vector wave equations contain two components, which 
are either Ex and Ey or Hx and Hy. This reduces the solution of the guided mode 
problem to two simultaneous equations for Ex and Ey or Hx and Hy. Using the 
vector wave equation for H, a similar wave equation can be derived similar to that 
for E discussed here. Formulating the problem in this way results in identical eigen- 
values as for the E expansion; however, the Et field results are inaccurate because 
the discontinuity of the E field is represented by a continuous H field [16]. 

To obtain a solution to the scalar wave equation. Ex is expressed as an 
infinite series in terms of a complete set of orthonormal functions. 


where 


oo oo 


y)) 

/1=1 v=l 


(4.111) 


= 


^Lxhy 


sin {(Ty,x) sin {p^y ) , 


(4.112) 
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is a simple and complete orthonormal set of sine functions, where 




/ITT 

r ; 

Jx 


and 


I/7T 

Pv — 7~ ) 


(4.113) 


(4.114) 


and fj, and u are integers. The sine functionj^'^i/, is orthogonal over the finite base 
domain of (0 < a; < Lx) and {0 <y < Ly), such that 

Lx Ly 

(4.115) 


//- 


0 0 


where 6 is the Kronecker delta function. For this method, the field is described only 
over this finite base domain; thus, proper choice of the domain limits is critical for 
this method to produce accurate results. The values for and Ly are chosen such 
that the field on those boundaries is zero. For well-guided modes these boundaries 
are small and this method performs well. However, for waves at or near cutoff, the 
field is not well-confined and the finite domain does not apply. 

For the Fourier method, equation (4.111) is inserted into the scalar wave 
equation, equation (4.60). Multiplying the left side by ^^j,(o:,y) and integrating over 
the finite base domain, produces a matrix equation expressed as 


mx \ 

I V J 

/i=li/=l \^o / 


(4.116) 


where 


Lx Ly 




LxLy 


r ty 9 


0 0 L 


dx dy, (4.117) 
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which can be rewritten as 


Lx Ly 


M, 




Lx Li 


J J n^{x,y)'<S^y{x,y)'^f,iyi{x,y) dx dy 


0 0 




Lx Lit 






(4.118) 

The index of refraction profile in equation (4.118) is approximated by rectangles of 
constant index, to obtain a closed form solutions of the integral. Thus, the rib was 
broken into rectangular regions of constant refractive index for the calculations [35]. 

Equation (4.116) can be rewritten in terms of a matrix C, with eigenvector 
^ . Increasing the number of terms, by increasing the limit 
values irix and rriy, increases the accuracy of the results. For many cases, a series 
of 10 sine functions in both directions is sufficient to produce accurate results [16]. 
In this study, 30 was used for rUx and rriy values, due to the large window size. 
This meant that the order of the matrix was 900, hence there were 900 eigenvalues. 
These eigenvalues represent both the guided and unguided modes, however, only the 
fundamental guided mode was of interest. 

Using the procedure from Marcuse [16], two FORTRAN programs [86] were 
developed to determine the propagation characteristics of the rib and fiber. The Lx 
and Ly chosen values, 15 fim and 10 fim, respectively, were large compared to the 
fiber and rib guide geometry. This choice was necessary to maintain the zero field 
condition on the domain boundary while the rib guide was shifted with respect to 
the fiber, in order to determine the coupling efficiency at various transverse offsets. 

Fiber Field - the Fourier method 

The Fourier method was applied to the elliptical-core fiber using the develop- 
ment just described to calculate both the field profile and the propagation constant. 


X and eigenvalue 
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The fiber dimension was closely approximated as a single rectangle with the inte- 
gration performed over the elliptical-core boundary. The same Lx and Ly domain 
boundaries as well as expansion terms used for the rib were used for the elliptical-core 
fiber. 

The Fourier Method - Validity 

The accuracy of the Fourier method and appropriateness of choices for the 
evaluation boundary and number expansion terms were studied. To determine the 
validity of the evaluation box boundary and number of expansion terms chosen, an 
analysis of the Fourier method was conducted to compare the E^(a:, y) values at a 
particular point for the rib guide. 

The Lx and Ly boundary values were varied from a 7 fim by 2 fxm box to a 
15 by 10 fim box increasing both dimensions in 1 fxm increments. Effects on E^(x,y) 
at (^- 0.05, f 

of the boxes produce similar E^ values with a maximum difference between points of 
5.99 X lO”"^, hence a 15 by 10 fim box is more than adequate for the structures being 
analyzed in this dissertation. 

For the rib, the mx and ruy values were varied from 5 to 30 in increments of 
5. Results shown in Figure 4.19 indicate that using 15 or more terms produces the 
same E^(x, y) value, hence the 30 expansion terms used are more than sufficient. 

Another analysis was performed to determine whether the electric field dis- 
continuity at the dielectric interfaces was evident. Field results taken at a; = 0 and 
^ ^ slices do not show the discontinuity, as illustrated in Figures 4.20 and 4.21. 

This is most likely an artifact of approximating the vector wave equation by the 


^ are shown in Figure 4.18. From this figure, it appears that all 
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Figure 4.20: Rib waveguide field magnitude 
at re = 0 vs. y calculated using the Fourier 
method. 

scalar-wave equation, but may also be due to the number of expansion terms used or 
the small dielectric discontinuity. 

As a proof of the accuracy of the Fourier method, it was applied to two 
different rib waveguide geometries found in the literature and field results were com- 
pared. The first rib. Rib 1, consisted of a 3 layer guide [60] as illustrated in Figure 
4.22. The calculated effective index was 3.437397 and 3.436799 using the EIM and 
Fourier methods, respectively. This rib was narrow and wide, hence the EIM is not 
expected to produce accurate results. Using the EIM and finite difference method, 
field profiles for the rib were given in the literature. The vast difference between the 
profile calculated using the EIM and finite difference illustrate how approximation 
techniques must be used with caution, because the EIM is not applicable to all rib 
guides. Contour plots of the literature results are shown in Figures 4.23 and 4.24. 
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Figure 4.21: Rib waveguide field magnitude at 
X = ^ \s. y calculated using Fourier method. 

Field results using the Fourier method, illustrated in Figure 4.25 closely resemble the 
finite element. 

Rib 2 differs slightly from a standard rib in that the index of the regions 
external to the rib are not of uniform index as illustrated in Fig. 4.26. The operating 
wavelength was 1.335 y.m and the rib width was 8 fim. The Fourier method predicted 
an effected index of 1.57058 which compares well to 1.57081 calculated using the 
domain integral equation analysis. Field intensity plots appear similar as shown in 
Figures 4.27 and 4.28. 

The propagation characteristics of a circular-core fiber was determined using 
the three solution techniques and results were compared to the exact solution given 
in the literature. The fiber and the rectangles used to approximate it for the Fourier 
method are shown in Figure 4.29. Effective index results are given in Table 4.4 and 
field results in Figure 4.30. The exact solution field result is shown in Figure 4.31. 
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Figure 4.26: Geometry of rib 2 from Kolk paper. 

Comparing Figures 4.30 and 4.31, it was evident that the Fourier field result was 
nearly identical to the exact solution. 

4.3 Coupling Efficiency Calculation 

There are several causes of interconnect loss between a fiber and a waveg- 
uide [50], [55], [57]. These include mode field mismatch, transverse offset {x,y offset), 
longitudinal separation {z offset), angular offset, and reflection. The loss mechanisms 
addressed in this study were the transverse offset and mode field mismatch. Mode 
field mismatch contributes significantly to interconnect loss [37], [57], hence modeling 
is a necessary tool to design a waveguide whose field profile closely matches that of 
the input fiber. Physically, mode field mismatch is directly related to the difference 
between the field profiles of the two components. Since misalignment in the x and 
y directions is inevitable in practice, the theoretical coupling efficiency results for x 
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Electric Field Magnitude vs. Rodiol d 



Rodial Distance (microns) 
Circular-core fiber 


Figure 4.30; Field amplitude versus radial dis- 
tance for circular-core optical fiber. 




Electric Field Amplitude vs. Radiol Distance Calculated Using 

Exact Solution 



Figure 4.31: Exact field solution for circular-core fiber. 
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and y offsets are of interest. Because the index of refraction difference between the 
fiber and waveguide was small, FVesnel reflection loss was neglected [57]. 

Once the electric fields of the rib guide and fiber are calculated, the efficiency 
of coupling between the two components can be determined. The coupling efficiency, 
rj, is calculated using the overlap integral [87], 

V 

where the integrals are taken over all space and E/{x, y) and Eg{x,y) are the electric 
field profiles of the fiber and rib waveguide, respectively. For the EIM and Marcatili’s 
method, this integration is performed directly. Recall for the Fourier method, both 
the fiber and rib fields are calculated using the same orthonormal basis functions, 
boundaries, and number of expansion (m and n). Thus, the coupling efficiency is 
determined as the integral of the product of the fiber and rib guide fields squared, 
which is 

(4.120) 

where and Bf^t, are the coefficients of the fiber and rib guide, defined by equation 
(4.116), which are normalized as 

= = (4-121) 

The coupling efficiency is related to coupling loss by the relation [60] 

Loss = —10 login 7? (4.122) 

The major axis of the elliptical-core fiber was oriented parallel to the x axis 
of the guide. The position of zero x and y offset is chosen to correspond to alignment 


j J Ef{x,y)E*{x,y) dx dy 

J J Efi^,y)Ef(x,y) dxdy J j Eg{x,y)E*{x,y) dx dy 
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of the maximum field positions of the fiber and rib [60]. After solving for the aligned 
case, various x offsets were introduced and a new coupling efficiency calculated for 
each at two temperatures. This procedure was then repeated for various y offsets 
while the x offset was zero. 

4.4 Concluding Remarks 

The field solution for a three-layer planar guide was discussed to provide a 
foundation for understanding the field solutions for the rib guide. Because there is no 
exact analytical solution for the rib propagation characteristics, three solution tech- 
niques were used. These were the EIM, Marcatili’s approximation, and the Fourier 
numerical method. Results from the three methods will be compared in Chapter V. 

The validity of the Fourier method was studied by comparing known solu- 
tions to those calculated using the Fourier method. Good agreement between the 
results was found, confirming the accuracy of the method and supporting its use as 
the baseline numerical method for comparison with the approximation techniques. 

The fiber does have an analytical solution which cannot be determined in 
closed form, however, it was determined that the three techniques used to solve the 
rib were also applicable to the fiber, and were expected to produce accurate field 
results. 


Table 4.4: Effective Index Results 


Method 

Effective Index 

EIM 

1.44931 

i^erturbation to Marcatili’s Technique 

1.44813 

the Fourier Method 

1.44817 




CHAPTER V 


RESULTS 

This chapter is divided into two main sections; room temperature results and 
300° C results. These two temperatures were chosen to represent a nominal and high 
temperature between which the guide is expected to operate when located near the 
inlet of an aircraft engine. Results from the coupling efficiency calculations performed 
at these two temperatures indicate that the device is expected to function properly 
in this environment. 

The field results for both the rib waveguide and elliptical-core fiber will 
be discussed. The field profiles and coupling efficiency results for x and y offsets 
will be shown in graphical format. Results for the effective index of both the rib 
guide and fiber as well as the maximum coupling efficiency calculated between the 
two components, using the three techniques, will be given in tabular format for 
comparison. 

5.1 Rib Waveguide Results (20°C) 

5.1.1 Field Results 

Theoretical field results using the EIM, Marcatili’s approximation and the 
Fourier method are illustrated using contour plots. For all the contour plots in this 
dissertation, each contour indicates a region of equal field magnitude in increments 
of 0.1, with the maximmn field magnitude being equal to 1.0. 
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The field profile calculated using the EIM is shown in Figure 5.32, with the 
rib cross-section superimposed. For this method, the solutions to Nr and Nn for the 
constituent guide were 1.5022 and 1.4938, respectively. These values were used to 
calculate the field as previously discussed. 

Using Marcatili’s method, the results for (3i and /?2 were 11.7134 and 11.3120 
1 ^, respectively. The contour plot of the electric field calculated using these variables 
is shown in Figure 5.33. Field results calculated using the Fourier numerical technique 
are shown in Figure 5.34. 



Figure 5.32: EIM: electric field contours of rib guide. 


5.1.2 Effective Index Results - Rib 

The results of the calculated effective index, Ug, for the rib waveguide using 
the three methods are given in Table 5.5. The result using Marcatili’s method was 
1.4963, which was improved by including the perturbation correction, 0.2186. 




METHOD n 


EIM 1.4994 


Perturbation (to Marcatili’s Method) | 1.4966 


Fourier 1.4993 
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5.1.3 Comparison 

Field results using the EIM were very similar to those from the Fourier 
method, yet required considerably less computation time. The difference between 
the fields calculated using the EIM and Fourier method are shown in Figure 5.35. 
Marcatili’s method assumed a tightly confined field, hence the results produced a 
field more tightly confined than the EIM and the Fourier method. The difference 
between the field calculated using Marcatili’s method and the Fourier method is 
illustrated in Figure 5.36. Results using the two approximation techniques were 
completed in a matter of minutes on a computer, however, using the Fourier method, 
the computation took a few hours on a Pentium computer. 

Effective-index results using the EIM and the Fourier method were nearly 
identical, with a difference of 0.007%; whereas between the perturbation method and 
Fourier method it was 0.18%. As evidenced here, the EIM is known to consistently 
predict a higher effective index than the actual value while the perturbation correction 
to Marcatili’s method predicts a lower effective index. Hence, the actual effective 
index value is bounded above by the EIM result and below by Marcatili’s result. 

5.2 Fiber Results (20°C) 

5.2.1 Field Results 

Modeling the elliptical fiber as a rectangle, the field results using the EIM 
are shown in Figure 5.37 with the fiber elliptical-core cross-section superimposed. For 
this technique the calculated Ni value was 1.4703 and 1.45 was used for Nn. Utilizing 
Marcatili’s method, the fiber field profile was determined using the calculated fix and 
/?2 values, of 11.1892 and 11.1300 respectively. A contour plot of the electric field 




Figure 5.35: Difference between Fourier field 
result and EIM field result for the rib waveg- 
uide. 



Figure 5.36: Difference between Fourier field 
result and Marcatili’s field result for the rib 
waveguide. 
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is shown in Figure 5.38. Field results using the Fourier method for the elliptical-core 
fiber are shown in Figure 5.39. Analogous to the procedure used for the rib, a 15 /im 
by 10 fim evaluation domain was used. 



Figure 5.37: EIM: electric field contours of elliptical-core fiber. 


5.2.2 Effective Index Results - Fiber 

The effective index of the elliptical-core fiber was calculated using the three 
solution methods with results shown in Table 5.6. The perturbation correction was 
0.3678, producing a corrected effective index of 1.4651. Using Marcatili’s technique, 
the calculated fiber birefringence was 2.5 X lO''^ [41] using equation (3.10), where 
and [3y denote the propagation constants for the TE and TM modes, respectively. 

This birefringence compares favorably with the manufacturer’s nominal value of 1.5 x 
lO--^ [78]. 




Fourier Method 
Elliplicol-Core Fiber (20 C) 



Figure 5.39: Fourier method: electric field profile of elliptical-core fiber. 


Table 5.6: Fiber Effective Index Results (20°C) 


METHOD 


rie 

EIM 


1.4653 

Perturbation (to Marcatili’s Method) 


1.4651 

Fourier 


1.4652 
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5.2.3 Comparison 

Both the field and effective-index results from all three techniques agreed 
well for the fiber. Differences between the electric field profile calculated using the 
EIM and Marcatili’s method with respect to the Fourier method are shown in Figures 
5.40 and 5.41. From these graphs, it appears that the field calculated using the EIM 
more closely matches that using the Fourier method than the field calculated using 
Marcatili’s method. 

The difference between the effective-index calculated using the EIM and the 
Fourier method was 0.0068% which was the same as the percent difference between 
the perturbation correction to Marcatili’s method and the Fourier method. The EIM 
overapproximated the effective index while the perturbation correction to Marcatili’s 
method underestimated it. In this case, the two approximation methods produced 
an effective index that was nearly identical to that of the Fourier method. 

5.3 Coupling Efficiency Results 

The zero x and y offset position was defined as the alignment of the maximum 
field positions of the fiber and rib such that the maximum coupling efficiency occurs 
at this position. An offset in the x direction was then introduced which ranged 
from -3 fim to -f3 fxm in 0.10 fxm increments. The coupling efficiency between the 
fiber and rib guide, calculated using the three different methods, is shown in Figure 
5.42 as a function of the x offset. It was evident from this graph that the EIM 
and the Fourier method predicted similar results. Marcatili’s method predicted the 
fields of both components more tightly confined in the x direction, which produced 
a coupling efficiency more sensitive to an x offset. For instance, given a 1 fxm x 
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Figure 5.40: Difference between Fourier field 
result and EIM’s field result for the elliptical- 
core fiber. 

offset, the predicted coupling efficiency using the Fourier method was 55%, whereas 
for Marcatili’s method it was 42%. 

Similarly for the y direction, keeping the x alignment, the fiber position 
was varied from -3 fjtm to +3 fim from the position of alignment. The theoretical 
coupling efficiency results versus y offset are shown in Figure 5.43. For a 1 fim y 
offset, the predicted coupling efficiency was 18.0% and 18.5% for the Fourier and 
Marcatili methods, respectively. The coupling efficiency versus y offset was nearly 
identical for the three methods. 

From the graphs in Figures 5.42 and 5.43, it was apparent that alignment 
was more critical in the y direction. This was expected since the rib field calculated 
was more tightly confined in the y direction than the x. Using the Fourier method, a 
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Figure 5.41: Difference between Fourier field 
result and Marcatili’s field result for the 
elliptical-core fiber. 



COUPLING [FFICENCY(%) 
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COUPLING EFFICIENCY vs. x OFFSET 



Figure 5.42: Calculated coupling efficiency vs. x ofTset. 


COUPLING EFFICIENCY vs. y OFFSET 



Figure 5.43: Calculated coupling efficiency vs. y offset. 
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1 iim X offset had a predicted coupling efficiency of 55% whereas a 1 y offset had 
18.0%. Offsets in the y direction produced coupling efficiency results that were nearly 
symmetric in y. Only a slight asymmetry in the coupling efficiency results versus y 
offset was determined, which was so slight that it could not be visually detected from 
the graphical results. 

The maximum theoretical coupling efficiency results are shown in Table 5.7 
to four significant figures. This was done to illustrate the minimal changes in the 
predicted coupling efficiency with temperature. Results using Marcatili’s method 
appeared to be the least accurate due to the method’s assumption of a tightly bound 
field. The results using the EIM and the Fourier method were similar, leading to the 
deduction that the EIM is the more accurate approximation technique for this type 
of rib structure. 

The zero offset position, in which the maximum field position of the fiber 
and guide are aligned does not necessarily produce the maximum coupling efficiency 
[88]. However, due to the symmetry of the rib structure and fiber being studied, the 
zero offset position was the maximum coupling efficiency position. This was evident 
in the symmetry of the coupling efficiency graphs versus offset. 

5.4 Rib Waveguide Results - 300°C 

5.4.1 Physical Changes Due to Temperature 

Compared to the 20°C7 data, the rib and fiber parameters changed slightly 
when calculated at 300° C. The nominal values for all the parameters were taken at 
room temperature. The physical features of the rib were determined at 300°C using 
the relations from Table 3.2, with the subsequent results given in Table 5.8. 
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Table 5.7: Maximum Coupling Efficiency (20° C7) 


METHOD 

MAXIMUM COUPLING EFFICIENCY 

EIM 

85.0348% 

Marcatili 

88.1608% 

Fourier 

81.9948% 


Table 5.8: Rib Guide Parameter Changes with Temperature 


PARAMETER 

Equation to Determine Change 

RESULT 

W{300°C) 

W{20°C) + (6.295 X 10"®^ x 280°C) 

2.701761 fxm 

^oi;er/a2/er (300 C7) 

^aueWaj/er(20°C7) + ^ — 9.119 X 10”®^^ X 280°C7^ 

0.508975 iim 

^tinetc/ied (300 C) 

tunetched{20’^C) + (-9.245 X lO"*^ X 280°C7) 

0.515974 fim 

tetched(S00°C) 

tetched(20°C) + (-7.632 X 10-8^ X 280°C) 

0.425979 urn 

'^overlayer (300° C7) 

n^erlayer{20°C) + (1.20345 X lO'^^ X 280°C) 

1.45637 

nni(300°C) 

tri6(20°C) + (1.1865 X 10-®^ X 280°C) 

1.55632 

'^isolation (300° C7) 

nisoiation{20°C) + (1.2035 x 10“®^ x 280°C) 

1.45637 
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The isolation layer thickness was not included in this table because that 
layer was modeled as infinite, hence any minute change in thickness would not alter 
the results. From the results shown in Table 5.8, it was apparent that the substrate 
stretched the film causing the rib width to increase while its height decreased. Also, 
the index of all three layers increased with the increase in temperature. The temper- 
ature induced change in the index of refraction was two orders of magnitude greater 
than the change due to the photoelastic effect, thus being the dominant effect. 

5.4.2 Field Results - Rib 

For the EIM, the rib field calculated at 300°C is shown in Figure 5.44. Field 
results for the rib waveguide using the Fourier and Marcatili’s methods are shown in 
Figures 5.45 and 5.46. Comparing these graphical results, the fields calculated using 
the EIM and the Fourier method were remarkably similar. These field contours were 
nearly identical to those calculated using 20°C parameters. 



Figure 5.44: EIM: electric field contours of rib at 300°C. 





Figure 5.46: Fourier method: electric field contours of rib at 300°C. 
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5.4.3 Effective-Index Results - Rib 

The effective indexes predicted by the three methods for the rib are given in 
Table 5.9. For the EIM, the values calculated for Ni and Nn were 1.5053 and 1.4968, 
respectively. The propagation constants, fix and calculated using Marcatili’s 
method were 11.7388 and 11.3876 respectively. The perturbation correction to 
the propagation constant calculated using Marcatili’s technique was 0.2182. 

5.4.4 Comparison - Rib 

Field results from the three methods at 300°C were nearly identical to the 
field results at room temperature. There was minimal change in the field results 
calculated at the two temperatures. 

The difference between the effective-index results calculated using the EIM 
and the Fourier method was 0.0003% whereas between Marcatili’s method and the 
Fourier method it was 0.15691%. This difference was less than that observed between 
the room temperature results. The effective index change between room temperature 
and 300°C was 0.00308 using the Fourier method. 

5.5 Fiber Results - SOO^C 

5.5.1 Physical Changes Due to Temperature Increase 

Using the values for parameter change per given in Table 3.3, changes in 
the indexes and core dimensions were determined and results shown in Table 5.10. 
Using these calculated parameter values, the field and propagation constant were 
determined using the three methods at this higher temperature. Compared to the 
room temperature values, the diameter of the fiber’s major axis increased by 0.004296 
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Table 5.9: Rib Effective Index Results (300°C) 


METHOD 


rig 

EIM 


1.50239 

Perturbation (to Marcatili’s Method) 


1.50003 

Fourier 


1.50238 


Table 5.10: Fiber Parameter Change With Temperature 


PARAMETER 

Temperature Change 

PARAMETER (300° C) 

a 

a(20°C) + (7.15 X 10-7°C x 280°C) 

1.302148 iim 

b 

b{20°C) + (3.58 X 10-7°C x 280°C) 

0.651074 fim 

'^core 

ncore(20°C) + (1.28 X lO-^* x 280°C) 

1.4868 

^cladding 

nc/<ufd.n,(20°C) + (1.28 X 10“^ x 280°C) 

1.4528 


Table 5.11: Fiber Effective Index Results (300°C) 


METHOD 


Tig 

EIM 


1.4681 

Perturbation (to Marcatili’s Method) 


1.4679 

Fourier 


1.4680 
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lim and the minor axis increased by 0,002148 fim. The indexes of both the fiber core 
and cladding increased by 0.0028. 

5.5.2 Field Results 

Field results for the fiber at 300° C calculated using the three methods are 
illustrated in Figures 5.47 - 5.49. Graphical results appear very similar to those at 
room temperature. For Marcatili’s method, results for /?i and P 2 were 11.2105 and 
11.1514 respectively. 



Figure 5.47: EIM: electric field contours of fiber at 300°C. 


5,5.3 Effective Index Results - Fiber 

The effective indexes predicted by the three methods for the fiber are given in 
Table 5.11. For the EIM, the calculated values of Ni and Nn were 1.4731 and 1.4528, 
respectively. The perturbation correction to the propagation constant determined 
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Morcatili’s Method 
ElMicor-Core Hber (300 C) 



Figure 5.48: Marcatili’s method: electric field contours of fiber at 300° C. 

using Marcatili’s technique was 0.3672. At 300°C, the difference between the results 
from the EIM method and the Fourier method was 0.0044%. 

5.5.4 Comparison of Fiber Results 

Results from the Fourier method were expected to be the most accurate. 
However, the computation time required by this method was orders of magnitude 
greater than that required by the approximation techniques. From the two approxi- 
mation methods, the field predicted by Marcatili’s method was expected to be more 
accurate than the EIM. However, graphical field results indicate the opposite, as both 
approximation methods had a 0.0068% difference in the effective index compared to 
the Fourier results. Compared to the room temperature result using the Fourier 
method, the effective index of the fiber increased by 0.0028. 
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Figure 5.49: Fourier method: electric field contours of fiber at 300°C. 
5.6 Coupling Efficiency Results 

The maximum field positions of both components were aligned for the zero 
offset position at SOO^C. Due to the symmetry of both components, the maximum 
coupling eflficiency occurred at this perfect alignment position. Graphs of the coupling 
efficiency versus x and y offsets are given in Figures 5.50 and 5.51. 

For the x offsets, Marcatili’s method and the Fourier method predicted signif- 
icantly different results. For instance, given a 1.0 (im x offset, the coupling efficiency 
was predicted to be 42% and 55% for Marcatili’s and the Fourier method, respectively. 
These are the same coupling efficiency values determined at room temperature. 

For the y offsets, the three methods calculated nearly identical coupling 
efficiencies. This was due to all three methods predicting similar y field confinement. 
For a 1.0 (xm y offset, the predicted coupling efficiency was 17.9% and 18.4% for 
the Fourier and Marcatili methods, respectively. This result was slightly less than 
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the 18.0% and 18.5% coupling efficiency predicted at room temperature for the same 
offset. 


COUPLING EFFICIENCY vs. x OFFSET 



Figure 5.50: Calculated coupling efficiency vs. x offset at SOO®^. 

Results from the three techniques are given in tabular form in Table 5.12. 
From this theoretical study, it appears that the coupling efficiency decreases slightly 
with temperature, by approximately 0.02% using the Fourier method. Results from 
Marcatili’s method predicted a slight decrease in the maximum coupling efficiency, 
0.01%, while the EIM predicted an increase of 0.2%. These predicted changes in the 
coupling efficiency with temperature are negligible and would be difficult to detect 
in practice. 

Table 5.12: Maximum Coupling Efficiency (300°C) 


METHOD 

MAXIMUM COUPLING EFFICIENCY 

Effective Index 

85.1815% 

Marcatili 

88.1462% 

Fourier 

81.9774% 
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5.7 Coupling Efficiency Optimization 

The optical fiber is commercially available, thus, its parameters cannot be 
altered. The waveguide can be altered slightly, within physical constraints, to maxi- 
mize the fiber to rib guide coupling. A preliminary study was performed to determine 
if the coupling efficiency could be improved. Each parameter was varied to deduce 
its affect on the coupling efficiency. It was determined that the coupling efficiency in- 
creased with decreased rib width and also with increased etch depth. These modified 
rib waveguide dimensions more closely match the fiber dimensions. Combining the 
two parameter changes produced the largest increase in the coupling efficiency, 6%, 
for a maximum coupling efficiency of 88%. Using these new parameters, the coupling 
efficiency versus x and y offsets was calculated. Results of this calculation are shown 
in Figures 5.52 and 5.53. Offset affects on the coupling efficiency indicate that the 
device is more lenient in the x offset tolerance while more sensitive to y offsets. 

5.8 Summary 

Results indicate that the EIM was surprisingly accurate for both the fiber 
and rib waveguide, given its lack of enforcement of the field boundary conditions. 
The EIM is known to consistently predict a higher ng, which was observed in this 
study. Temperature induced variations of the rib and fiber indexes and dimensions 
were much less than expected. Due to the different coefficients of thermal expansion 
for the rib layers, the rib was stretched by the substrate expansion which slightly 
decreased its height. The indexes of both the fiber and rib guide increased slightly 
with the temperature increase. The field profiles of both components at the two 
temperatures were nearly identical. While results from the Fourier method were more 
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COUPLING EFFICIENCY vs. y OFFSET 



Figure 5.51: Calculated coupling efficiency vs. y offset at 300°C. 


COUPLING EFFICIENCY vs. x OFFSET 
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Figure 5.52: Optimum rib waveguide coupling 
efficiency vs. x offset from the preliminary 
study. 
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COUPLING EFFICIENCY vs. Y OFFSET 



Figure 5.53: Optimum rib waveguide coupling 
efficiency vs. y offset from the preliminary 
study. 

accurate, the amount of time required to determine the coupling efficiency for x and 
y offsets was several orders of magnitude greater than that for the approximation 
methods. 

These results indicate that the device is likely to perform properly in an envi- 
ronment with a temperature fluctuating between the two temperatures investigated. 
The calculated maximum coupling efficiency was higher than expected. Changes in 
the coupling efficiency with temperature were negligible, and the parameter changes 
with temperature were lower than expected. 




CHAPTER VI 


CONCLUSION 


6.1 Summary 

The electric field profiles of an elliptical-core fiber and a SiON rib waveguide 
were determined using three methods: two approximation and one numerical. These 
calculations were performed for room temperature and 300°C parameters. From 
the field results, the maximum coupling efficiency between the two components was 
calculated from the overlap integral and found to be higher than expected. Prior to 
this work, the coupling efficiency had been calculated between a single-mode circular- 
core fiber and the rib waveguide under study. Results from that work indicated an 
expected coupling efficiency of only 45% [75]. Therefore, choosing an elliptical-core 
fiber improved the device coupling efficiency by nearly a factor of 2. 

These results indicate that the device is likely to perform in an environment 
whose temperature ranges from 20°C - 300°C. Changes in the physical parameters 
with temperature were lower than expected over the 280° C range. 

As an offshoot to the problem under study, a comparison of the three solution 
techniques was performed. The Fourier method was expected to produce the most 
accurate results for both the fiber and rib guide. Therefore, compared to the Fourier 
method, it was determined that the EIM produced more accurate results for the rib 
guide and elliptical-core fiber. For an additional comparison, Marcatili’s method was 
to determine the fiber field and the EIM for the rib. The predicted coupling efficiency 
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was 81% at room temperature [89], which compares favorably with the Fourier pre- 
diction of 82%. Calculations using the Fourier method are computationally intensive 
while results using the approximation techniques are achieved in minutes, making 
the approximation techniques useful for preliminary studies. 

The theoretical coupling efficiency between an elliptical-core fiber and a rib 
guide was determined for a range of transverse offsets. Three techniques were used 
to calculate the theoretical electric field profile of the rib waveguide: the EIM, Mar- 
catili’s technique, and the Fourier method. For this channel structure, the EIM and 
the Fourier method produced similar results. Results from Marcatili’s method were 
not as accurate for the particular rib guide being studied. Similarly, for the fiber, the 
EIM was determined to be more accurate than Marcatili’s approximation technique. 

The coupling efficiency results determined using the three methods compared 
well. Calculations of this kind are essential to optimize the design of a rib guide in 
order to provide an acceptably low coupling loss. 

6.2 Significance of Work 

The fiber and rib waveguide were determined to be compatible for intercon- 
nection over the 20°C - 300 °C temperature range. Minimal changes were determined 
over this temperature range for the field profiles of the components and the associ- 
ated coupling efficiency. The rib waveguide under study could be implemented as 
the input stage for a variety of sensors fabricated using SiON as the guiding layer. 
These sensors would be expected to function in the temperature environment of the 
inlet of an aircraft engine. Coupling a single-mode fiber to a single-mode device is 
more sensitive to offsets than coupling a multimode fiber to the same device, due 
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to the smaller diameter of the single-mode fiber. However, these results indicate 
that greater than 50% coupling efficiency is predicted for x offsets within ± 1 {xm 
and y offsets within ± 0.5 fim. Offsets in the y direction were determined more 
critical to the coupling efficiency than in the x direction, which must be considered 
when designing the coupling configuration. These results are encouraging for future 
experimental research on interconnection for this fiber and rib waveguide. 

The analysis also provided a comparison between three techniques with re- 
sults from the Fourier method expected to be the most accurate. This demonstrated 
the surprising accuracy of the EIM for both the rib and fiber. Preliminary calcu- 
lations can be performed using these approximation techniques without sacrificing 
much accuracy, which will greatly decrease the amount of computation time required. 
Marcatili’s method was determined less accurate than the EIM for both the rib waveg- 
uide and fiber. Approximation techniques must be used with care because they are 
not applicable to every guide geometry. 

6.3 Future Work 

Plans to experimentally couple the light output of an elliptical fiber into the 
rib waveguide sensor at room temperature are under way. Differences between the 
theoretical and experimental results will be studied to determine the cause. Generally, 
the experimental coupling efficiency versus xovy offset matches the theoretical results 
multiplied by a transmission coefficient. This transmission coefficient is used to 
account for reflection loss, propagation loss, and output coupling efficiency [48]. 

In regards to the rib waveguide, a few researchers have reported a phe- 
nomenon where a rib which begins as single-mode will become multimode as the etch 
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depth increases, which is expected. However, as the etch depth further increases to 
a critical point, the device become single-mode again. This has been attributed to 
spurious modes caused by the finite window size. Modeling the rib and varying the 
etch depths to replicate this phenomenon would be of interest. 

Calculations to determine how the coupling efficiency changes with varia- 
tions in the gap between the fiber and guide would be useful for determining device 
operation. It would also be interesting to model the rib waveguide using four layers, 
instead of three, to determine any differences between the results. 

The design for coupling the fiber to the waveguide will use a v-groove configu- 
ration, hence thermal changes of the groove will affect the alignment. This alignment 
change with temperature was not considered in this research, however, during the 
design stage, these thermal effects will be considered and the subsequent coupling 
efficiency determined. To eliminate loss due to scattering at the fiber/rib interface, 
it will be necessary to polish the endface of the rib waveguide. 

The issue of bonding the fiber into the groove must also be addressed. Using 
epoxy as a bond is not a practical choice because epoxy degrades at high temperatures 
and has a different coefficient of thermal expansion than the fiber and guide, thus 
altering the alignment with temperature. A flip-chip approach may be applicable to 


this situation. 
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